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ABSTRACT 


The metiod of turbine performance prediction developed by 
Vavra and Eckert has been refined in this analysis to realize more 
of the potential of the three-dimensional calculating methods. 

Mach number and rotor tip clearance effects on blade outlet angles 
and loss coefficients have been localized rather than averaged over 
the blede height. An approximation for streamline curvature has 
been used. 

Performance curves were determined for two single stage axial- 
flow turbines located at the Propulsion Laboratory of the Naval 
Postgraduate School. Test results were available for one of the 
turbines. Agreement between predicted and experimental performance 


values was generally within 3 per cent, 
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Thesis by Robert G. Harrison entitled: "An Analysis of Single-Stage Axial- 
Flow Turbine Performance Using Three-Dimensionak Calculating Methods" 
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1. Introduction 

Turbines form an important part of propulsion systems. To opti- 
mize a design it is necessary to know the performance at off-design 
conditions as well as the performance at the deSign point. Since the 
testing of a prototype is very costly and time consuming, it is of 
great advantage to be able to predict the performance characteristics 
by means of theoretical methods. The more advanced a system is the 
more important it becomes to improve the accuracy of these methods. 
Hence it is necessary to base these methods on the fundamental laws of 
fluid dynamics rather than on ru.e-of-thumb approximations. The latter 
basis is possible only when data on previous designs are available. 
For new and advanced configurations, it will be necessary to apply re- 
fined methods which enable the designer and systems engineer to pre- 
dict the effect of proposed design changes. 

The principal equations that describe the flow properties in 
turbomachines are well known. These same equations are used as a basis 
for ail proper ‘'three-dimensional”’ calculating methods that have been 
cited in the technical literature. The methods differ however in the 
manner in which these equations are manipulated and applied. 

This thesis is concerned with the refinement of the three- 
dimensional method of analysis developed by Vavra and Eckert. Based on 
the physical dimensions of particular test turbines that are available 
at the Turbo-Propulsion Laboratory of the Naval Postgraduate School, 
performance curves were determined for these machines. Concurrently 
with the performance analysis, experimental tests were conducted on 
one of these turbines so that actual experimental results could be 
used to julge the accuracy of the proposed theoretical performance 
evaluation. 

In addition to his publications and classroom lectures which con- 
Stituted the foundation for this analysis, Professor Vavra was very 
generous in providing guidance and counsel during the period of this 
work, For this I am greatly appreciative. I would also like to thank 
Lieutenants P, M, Commons and J. A. Messegee for making their experi- 


mental results available. 
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2. Basis for the Analysis 

The two conservation equations that must be satisfied to obtain 
a solution for the flow in a turbine are the equations of motion and 
continuity. These equations are satisfied at stations between blade 
rows. The method used is that given by vavra.* Vavra developed the 
equations of motion and continuity for absolute flows. The equa- 
tions in this form are readily useable for the position after the 
stator. Eckert later developed these equations for relative flows to 
be used for rotor eafentlartons.” Eckert's conversion allows relative 
flow quantities to be handled directly without conversion to an 
absolute system. Eckert's approach also avoids iteration procedures 
to determine the total enthalpy after the rotor. 

The assumptions made for the developement of the equations 
used in the performance analysis are: 

1. An infinite number of blades in each row so that downstream 
effects are not felt upstream. 

2. Axisymmetric flow at the stations where the equations of 
motion are solved, 

3. Adiabatic and steady flow so that the total enthalpy along 
any given streamline is constant through the stator and the relative 
total enthalpy is constant through the rotor, 

4. All entropy changes are assumed to occur in the blade 
channels that are located ahead of the stations where the equations 
of motion are satisfied. Hence at the calculating stations the 
flow is assumed to be isentropic along particular streamlines. 

With the above assumptions, the equations of motion for 


absolute and relative flows, respectively, are 


JH= Vx(vxV)+T 95 we 


GHe=Wx(¥xW+20)+T 9S (2) 





yavra, M. H., Aero-Thermodynamics and Flow in Turbomachines. 
New York, London: John Wiley and Sons, Inc.; 1960, Chapter 16. 


anckert, R. H., Performance Analysis and Initial Tests of a 
Transonic Turbine Test Rig (USNPGS Thesis, May 1966), pp. 149-155. 
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Relative tct2i enthalpy can be written as 





2 2 2 2 - 2 
iy tee mall ~ We _ Ve - fb , Ween Ve (3) 
ail licens 2. 2 2. a 2. 2 
Equivalent erthatpy a aefined as 
Ween _ Us 4 
H_e= nee a = Hr oe ee 2 (4) 
Similar to Hs the enuivalent enthalpy Hi. is constant along a 
streamline for the adopted assumptions. Equivalent enthalpy can 
also be written as ae que oa 
4 =h Wes h, ._—— — Rd! (5) 
E” "oaa 2 


The introduction a the © enthalpy allows this quantity to 
be used for the rotor in a manner analogous to the way total enthalpy 
is used for the stetor. 

For equations used in this analy: is, the subscripts refer to: 

0 - station ahead of stator 

1 - station ahead of rotor 

2 - station after rotor 

is - isentropic expansion from re 
s - isentropic expansion from Ree 
th - the theoretical value ; 

Figure 1 is a temperature-entropy diagram showing the thermo- 
dynaric process along a particular streamline for a single stage 
turbine, In general, the fluid properties will vary from stream- 
line to streamline, Thre methow in which the loss coefficients are 
appliea is also incicated in Fig. 1. The loss coefficients are 


defined ~, 


ven Ms 
—— For the Stator (6) 
ae 
2 
ws Wee - We. For the Rotor (7) 
R- Ween 


The coordinate system that will be used in the analysis is 
indicated in Fig. 2. This figure also shows the general layout of 
the tv; turbine to which the prediction performance analysis is 
applied, 

ign convertion for the various arglies that are needed in the 


analy:.iz { indicated in Fig. 3. 
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The modification of Eq. 2 into a form that can be used for the 
analysis is given in Appendix A, Section 1. The appropriate form of 
Eq. 1 can be obtained from the modification of Eq. 2 if the angular 
velocity w is set equal to zero, Other differences of the final 
equations derived from Eqs. 1 and 2 are listed in Appendix A, 
Section l. sis 1 can then be can 





ULE) ery _/LzeRe)*) st dhe 
a 2 x,|(-K2 rm ae) os Le a  |- = ae 
Oe ee “x, dH mir. Heos*a, eins | ee. 

e SIN X, te % Var d ; | ve Vn | d X, (8) 


Equation 2 becomes 


dhe) “) E Le eR/ar ae, ee ee “ 
“dX, = ~cos'£, (K2 Te. or). SS L2 ra allel oar e3 ge 


_ +n cashes on Be _ _—_— + ses Ae Tie [Gl ce oat (9) 
Wares a Wa Ye dX2 L Waa. Yo “Id Xo 
where: 
a or a for Eqs. 8 and 9 respectively 
View ts (subscript m refers to mean streamline) 
= 
=a 
K=5.0 streamline curvature factor 
Sy streamline displacements shown in Fig. 2 
AR 
= Lithlke 
—— 3 = (see Fig. 2) 
_S 
>= "n 
eee) 


The equation of continuity is used in its non-dimensional form 


by introducing a flow function d . Developement of $ is given 


in Appendix A, Section 2. In differential form the equations for 


absolute and relative flows can then be expressed by 


ville [Es ang= aay SCE) (ES) | ao 


‘2 


and 


a aahadenaseemmmmenandamiane ee 


if [wands BRIE ° 


—E 
The differential element of area dA is 


JA= Ezadr (12) 
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where: 

2 = number of blades 

a = blade exit opening (see Fig. 46) 

§ = are: restriction coefficient 

Since $ is valid for isentropic flow only, the restriction 
factor 5 must be introduced to correct the actual flow area to an 
effective area which accounts for the restrictions due to the 
boundary layers on both sides of the flow channel. 


The factor € can be expressed her 


x wt 


ee ee (13) 
k; = i eee ‘ 
In Eq. 13, ¢ * is the koundary layer displacement thickness and 
hx** is the so-called energy parameter defined as the energy thick- 
ness civided by the displacement thickness. The term = repre- 
sents the loss that is assumed to occur from the inlet to the 
throat of the blade channel, where ‘/ is the loss coefficient repre- 
senting all the losses across the row of blades. The profile loss 
coefficient wes used by petents to represent the loss prior to the 
blade throat. Percentage of the total loss due to profile losses 
will vary considerably depending on blade geometry, radial position, 
and the incidence of the flow on the leading edge of the blade, 
Since secondary flow and tip clearance effects result in losses in 
the blade channel, half the total loss coefficient provides a better 
average representation of the losses in the blade channel prior to 
the throat. Thre basis for the developement of = and H*** as used 
is given in Appendix A, Sections 3 and 4. 


By multiplying and dividing by a, 1. Eq. 12 is 


i ee (14) 


2 pte ar 


wee 


“Yavra, M, h,, Problems of Fluid Mechanics in Radial Turbo- 
machines (i:hode-Saint-Genese, Belgium: Von Karmdn Institute for 
Fluid Lynamics, 1965) VKI Course Note 55b, pp. G46-50. 


4 
nckert, op, cit., p. 44, 
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After integration Eqs. 10 and 11 become, respectively, 


Kr 
WATro R _ Qa 


Sy =e Zr. fi Ed ik | (ey 


ae 
The flowrate w can Be computed from the conditions ahead of the 


and 


Stator. Then a reference flowrate is defined by 


{vite [R | 
Gla a | oy. 


where W eee is in square inches. Continuity will be satisfied for 


the stator by 


Ar 
E =. fa E g x = Wee; (18) 


R) STATOR 


Similarly for the rotor 


antrnf tof 8 eT dX 
tc "Cae 
= ROTOR 


‘the ae of the leakage flow through the radial tip clearance 





= Wres (19) 





has not been accounted for in Eq. 19, 


The element of area between the blade tips and the shroud is 
th=arrdr (20) 


the flow through the tip clearance area is 


a Xs 
M| T+: ae : x 
= bs andy Bt | X &4g dx oe 
te me 
IE TP hr 
WEAR ENCE 
Since the tip clearance is relatively small, the values of 
f t P. , T, , and X, for the tip will be used in Eq. 21. 
’ ra te tr Zz 


With these assumptions Eq. 19 a) be expressed by 


“ee Te a ee ern "tEs(Be EEX ) et ee 


Wres 


= (ane Ym] Ro 
The assumptions used to arrive at Eq. 22 are obviously incorrect in 


C225 


two respects. First, the flow represented by $ is not 
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perpendicular to the tip clearance area, Second, the effective area 
represented by EL is larger than that which probabiy occurs 
because of the relatively large boundary layers that exist on the 
Shroud and blade tips. The exact behavior of the flow in the small 
region between the shroud and the blade tips is impossible to predict 
without further tests. However, it is felt that the flowrate through 
this space as represented in Eq. 22 is too large for the reasons just 
mentioned, Therefore a more accurate approximation of this flowrate 


will be obtained if the last term on the left side of Eq. 22 is 


rai by 2, yielding Xs 
“bef 2 gg ax vom f(SfE g3x) | | 
ky, @ Re, ae 
ve TIP ROTOR 
= iret. (23) 
(Lm Z Val ROTOR 


The tip clearance flow included in &q. 23 can be obtained by inte- 


Saesar 


nt f RET Ho FE (EVE SPY) be 
WR ag 


—- +e wkir (Re |The 
x “(= Tee. eX) TIP Veemall - the | | Tee ef) (24) 


[neroduetng ¢ this expression into Eq. 23 gives 








Ko 
pete 2 pga ee (Saye 2), J 
Ky Re \ Tee Ca = ame Yin Fee | Tee ~ “/riel ROTOR 
_ Wre 
Thin (25) 


3. Technique for Obtaining Solution 

With the equations of motion and continuity in the forms given 
by Eqs. 8, 9, 18 and 25, a method has been developed to analyze 
Single stage axial turbines, in particular, those available for test 
in the Turbine Test Rig of the Turbo-Propulsion Labora ory of the 
Naval Postgraduate School. The method of analysis predicts turbine 
performance for specified values of inlet total pressure, inlet 
total temperature, rotor speed and the ratio of total inlet to static 
discharge pressure P IP. This method is similar to that described 


c 
by Eckert.> However, °Eckert's analysis neglected some effects which 


"Ibid, Section 3. 
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have been accounted for in this developement. Some significant 
changes made in this method are listed below: 

1. Stator and rotor outlet angles for a particular radial 
location are computed using a calculated Mach number for that loca- 
tion rather than an assumed Mach number or the Mach number of the 
mean Streamline, 

2. Variation of loss coefficients due to changes in blade geo- 
metry in the radial direction is accounted for. 

3. The influence of rotor tip clearance on the rotor outlet 
angles and loss coefficients is concentrated near the tip of the 
blade rather than averaging these effects over the full blade height. 

4. Fourth order polynomials are used to better approximate 
the curves respresenting blade characteristics as a function of 
radius and the curves of rotor loss coefficients as a function of 
incidence, 

5, Streamline curvature effects have been accounted for in the 
solution of the equations of motion, 

In addition to the assumptions that were mentioned in Section 
2 for the developement of the particular form of the equation of 
motion, the conditions ahead of the stator are assumed to be uniform; 
that is, the total temperature, velocity, and entropy are assumed 
constant and the flow axial in direction. It is realized that com- 
pletely uniform conditions are difficult to obtain, but any other 
assumpt on would be extremely difficult to develop mathematically. 

Direct solution of the equations of motion is not possible since 
they are nonlinear in the dependent variable Y, Likewise, no direct 
method is possible to satisfy continuity. Solutions of these equa- 
tions must therefore be gained by making initial assumptions for the 
values of the axial velocities which must be improved by successive 
iterations until the equations are satisfied, To account for stream- 
line curvature and slope, a complete solution of the flow through 
the stator and rotor must first be made by neglecting the effects 
of curvature in order to determine streamline locations, Then the 
iteration to account for these effects may progress, These require- 
ments make the use of a high-speed computer a necessity. 
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This analysis has been programmed for the IBM 360 computer 
using FORTRAN IV, The program is described in Appendix C. The 
following paragraphs set forth the procedural steps of the program. 
The equations are listed in general form without referring to 
specific streamline locations. In the interest of clarity, however, 
some relationships will be written in forms similar to those used 
in the program. For example, P(2/2,=[ Tis @)/Te JF will repre- 
sent the isentropic relationship for the number 2 streamline. 

Five streamlines are utilized for the analysis with the 
number 1 streamline located at the hub and the number 5 streamline 
located at the tip as shown in Fig. 2. The number 3 streamline will 
be used as the mean streamline, and the radius of this streamline 
will be referred to as the mean radius. The radial locations of the 
Streamlines ahead of the stator will be such that the mass flowrate 
between adjacent streamlines is 25 per cent of the total flowrate. 
Positions for the streamlines after the stator and after the rotor 
are initially assumed. The locations of streamlines 2, 3, and 4 
then vary during the solution as necessary so that the percentage of 
the total flowrate between adjacent streamlines does not change. 
This continuity requirement will be called streamline continuity. 

Besides the radii, sufficient input information must be used 
to effectively reflect the physical characteristics of the stator 
and rotor blading. Some of the physical properties are introduced 
directly; such as, the number of stator blades, the number of rotor 
blades, and the rotor tip clearance. The other quantities used which 
reflect blade characteristics are throat opening dimensions for the 
blade channels, discharge angles, rotor blade inlet angles, loss 
coefficients, and stalling incidences for the rotor. 

Throat opening dimension ''a" is a function of radius. The 
best method for introducing this characteristic into the anaylsis 
is to enter the measured values of "a" together with the correspond- 
ing radii. Then, utilizing the method of least squares, a fourth 
order polynominal curve is fitted through these points. From the 
resulting polynomial, the value of "a'' for any radius required by 


Streamline continuity can be determined, 
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Discharge angles are nredicted by using a combination of the 
methods of Vavra and Rinlevse Outlet angles are first calculated 
using the formula which Vavra established from the experimental data 
of 2 awesld These angles are then corrected for tip clearance, blade 
cu. vature, and Mach number effects with the methods given by Ainley. 
Stator and rotor discharge angles are predicted at three radii; 
namely, the hub, mean radius, and tip. The method used for calcula- 
ting these angles is explained in Appendix B, Section 1, 

Values of stator gas outlet angles &, for the mean stream- 
line are determined for Mach numbers My of 0.5, 0.7, 0.75, 0.8, and 


1,0, These values are represented by two parabolic curves of the 


form 

= a+ oe mill (26) 
The first curve is used for Mach numbers M, from 0.5 to 0.75 and the 
second established interim values of Xn, for M, between 0.75 and 
1,0, From these curves the flow angle &, for the mean streamline 


can be found for any Mach number M,. 


The flow angle ©, is also a .zunction of radius Tr): There- 


fore the changes of &, ‘or the hub and tip with reference to the 
mean radius, called 4%, and Air respectively, must be used. 
The flow angle &, can then be determined for any zy by assuming 
a linear variation or M, between the hub and the mean radius and 
between the mean and the ti», With this assumption and using the 


Mach number M, in Eq. 26 corresponding to the number 2 streamline, 


1 
the flow angle &, for this streamline would be 


Bee” Wale) - ae 
(X ee) = (x | ro — eee O”% Ze 
ce “— 


et SO ere Ee 


Sa iniey, D. G, and Mathieson, G, C, R., A Method of Performance 
“atimation for Axial-Flow Turbines. Aeronautical Research Council, 
h & M No, 2974, 1957. pp. 34, 


"Heer, k,, Aerodynamic Design and Estimated Performance of a 


Two-Stage Curtis Turbine for the Liquid Oxygen Turbopump of the M-)] 
Enpine, NASA CR 54764 (AGC 8800-12), Nov. 19, 1965. p. 29. 
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The superscript ** is used with r, and qh, in Eq. 27 to indicate 


the radius initially assumed for the mean streamline and the 
computed for that radius. Equation 27 can then be used throughout 
the analysis even though the radial location of the mean streamline 
may change due to streamline continuity requirements. A similar 
approach is used to establish the flow angle &, at the rotor dis- 
charge, 

The rotor blade inlet angles & are measured from the manu- 
facturing drawings of the blade profiles. Using the values of G, 
for the hub, mean, and tip streamlines, a parabolic curve is deter- 


mined which gives Z& as a function of xX): ; 


Loss coefficients and stalling incidences are predicted by 
using the methods of Ainteome For the present method, the stalling 
incidence i. is defined as that at which the loss coefficient is 
twice the value of the minimum loss coefficient. Following 
Ainley's methods, loss coefficients are computed as e function of 
the ratio of flow incidence to stalling incidence £ . Loss 
coefficients are also a function of blade geometry eee Since 
the stator has zero incidence, its loss coefficient is a function 
of the radius only. . 

Stator and rotor loss coefficients and rotor stalling incidences 
are calculated for the hub, mean, and tip radial locations. Rotor 


loss coefficients b are determined for values of ranging 


i 
CS 
from -2.0 to 1.6. Curves of ise vs. ae are drawn for each of 
the three radial locations, and values of ce ,» along with the 
corresponding quantities oe » are used to determine fourth 

order polynomials which approximate these curves. A similar proce- 
dur is followed to obtain a fourth order polynomial representing 
stalling incidence i, as a function of radius Yr, . Sample calcula- 
tions for the prediction of stator and rotor loss coefficients and 


stalling incidences are contained in Appendix B, Section 2. 


orale op. cit., pp. 4-5. 


Variation of loss coefficients with radial location is accounted 
for by assuming a linear variation of these quantities between the 
hub and the mean radius and between the mean and the tip. To demon- 
strate the procedure followed in computing rotor loss coefficients, 
the following example is given. For a particular incidence, the 
first step in the determination of we for the number 2 streamline 
would be to calculate by by using the radius Y, of that streamline 
and the polynomial of the form ins i (r)). Loss coefficients for 
the hub and mean radius would then be computed using the resulting 
~ in the polynomiale, for these radii, of the form Y= % "Ga 
with the assumed linear variation ‘ (2) would ve 


(a= (i) + tie cal lke 
The reason for using the suncrscript ** on the mean streamline 
values is the same at previously mentioned in connection with 
Ea. 27, 

For the first approximatior,the Mach number M, ahead of the 
stator is assumed,and the static properties and flowrate at station 


© are found by 





- Wi|Rle 2 w 
gA* ae lede” Sree 


The reference flowrate ee will be used to check continuity at 


stations 1 and 2. 
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(28) 


(79) 


(30) 


(31) 


(32) 
(33) 


(34) 


(35) 


The next step is to determine axial velocities after the 
stator that satisfy the ecuation of motion. Total enthalpy after 
the stator is constant by assumption, and streamline curvature 
effects are neglected at this stage of the analysis. With these 


conditions, Eq. 8 simplifies to 


2 2 * 
din Yi.) dz, 2 2 ( C,Hcos*d, \98 
= 2TAN @, = —- = SIN WA +t [= ze ao 
OX, ‘AX, =X, / v Va; m OX, (36) 
This equation can be integrated to give 
X 
InY*= ( TdX + Inc” 
Xo (37) 


where xX, is arbitrary and /nc* is the constant of integration. 


Using fie boundary condition Y = 1.0 at X = 1.0, Oe can be 


ime on InC*= - {Tdx 
Ke 


y2=f1 IdX- - SI TIX = St IdX x 


act cn 38 ee be eoeenee in a form that can be utilized 


found by 


then 


in the feo ton Expansion by infinite ne yields 


X 
2 A sales be | 


The quantities contained in I, Eq. 36, must be evaluated i 
before proceeding with the solution. For the first approximation, 
assumptions are made for the values of Yy M,> and Va . The flow 
angles Q, are then calculated by using Eqs. 26 and "07, After the 
value of Q, has been determined for each streamline, ses is 
computed. This derivative and all others needed in the analysis are 
found by finite difference methods. 

Enthalpy is computed by 

W= op Te, (40) 


and the entropy term is found by the method of Vavra. 


~ 


"Vavra, M, H,, Aero-Thermodynamics and Flow in Turbomachines., 
New York, London: John Wiley and Sons, Inc., 1960. pp. 445-447. 
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- ie Var me 
ae he CH Cos? (41) 
YP Vain 
a C,H cos*a (i- 4%) 
The stator loss coefficients Se are determined by equations similar 
to Eq. 28. 


Each time a solution to Eq. 36 is found, the new values of yy, 
are then used for the next iteration. After five iterations, &, 


and cal 





are recalculated using the new stator exit Mach numbers. 


\ 
The Mach number for each streamline is found by 














Va, =f it VA (42) 
V = Vay (42a) 
| Cos, 
Ni 
Ty, Tee 7 agscp (43) 
Vi 
= 44 
Mi Fegrt, als om 
With the new values for X, and — , five more iterations are 


made to determine the corrected valitied of Ya: 


The quantities represented by Eqs. 42-44 are recomputed after 
satisfying the equation of motion and additional quantities deter- 





mined by 
Vu, = Va, TAN & (45) 
ie | 
ese 7-5, (46) 
Th ee 


The flowrate through the stator is computed next and compared 
with that required to satisfy continuity. Reference flowrate between 


the hub and each streamline is found from 


XK X 
SUM = anf fn & E bdk = and rmf WI dX (48) 
in Ay 


where 


ssl et” ‘ 
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The "a" and ‘ae in Eq. 48 are found from the polynomial which 
represents throat opening as a function of radius, 

Before ug. 46 is solved, the pressure ratio is compared with 
the critical pressure ratio, if the critical pressure ratio has 
been exceeded, the flow is choked at that radial location and the 
critical pressure ratio is used in calculating p and = P 

The fraction of the total flowrate passing between the hub and 


each streamline is computed by 


es (50) 
WI 4x 
Total flowrate, as found by the denominator of Eq. 50 multiplied 
by uy, , is then compared with the reference flowrate. 
Overall continuity is satisfied if the difference is less than 
0.0002. 
If required and computed flowrates are not within this toler- 
ance, the axial velocity for the mean streamline is adjusted by 
Vlre:.7 Weempute 

Vaeeeen aise se) “ees Soo 

ree is the required reference flowrate divided by @,2Y,,, and 


, is the denominator of Eq. 50. 
computed 


Solutions to the equation of motion and continuity are 
successively found until overall continuity is satisfied. 

The fractions of the total flowrate determined for each 
Streamline by Lq. 50 are then compared with the corresponding flow- 
rate fractions ahead «cf the stator, Streamline continuity is satis- 
fied if agreement is within 0.002 for each streamline, If stream- 
line continuity has not been satisfied, the streamlines in error are 


adjusted by 


: c dX 
x We ny = boa +| Ween.” We oo Pied (52) 
Equation 52 applies only to stieamlines 2, 3, and 4 since We = 0 at 


the hub and W, = 1.0 at the tip. 


f 


If streamline positions have been adjusted, new streamline 


radii are found by 


fe 


Voty = Anew Y* ovo ; (53) 


and new values of X for each streamline are obtained by 


Vrey: 
X = ae (54) 
™ NEW 
With the new streamline radii the equation of motion and overall 


continuity must be satisfied again. 
The following values are determined after streamline continuity 


is satisfied: 


_ REM 
Ui Goyaay % se) 
i> 
U.= U, (56) 
Wis = Vu, ~ U, (37) 
gue =| (Nona 
[2 = TAN Gre (58) 
Vay 
Wi* 2o88, (59) 
2. 2 a 
Ww U2 7 
Ve es £9%Cp * Bqacp (60) 
Rem lbs CP (61) 
ce 7 1-1 
ee ie ee ) (62) 


Rotor blade inlet angles are determined for each streamline 
by means of the parabola which establishes veh as a function of Xy- 


Incidence is found by 


= Janes (63) 
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Stalling incidences and rotor loss coefficients are then deter- 
mined using the procdure described in connection with Eq. 28. If 


the incidence ratio is less than <-2.0 or greater than 1.6, 





Ss 
the value for —— is set equal to -2.0 or 1.6, respectively. 
( 


S 
One of the quantities necessary for the equation of motion after 
ds2 


d 


the rotor is This quantity is separated into two parts, 





X 
: dar dsei (64) 
AY, fi ¥ cal Xo 











| * 

C1819 represents the entropy gradient due to changes of entropy , 
cS 2) 

ad Ko 
represents the entropy gradient due to the different entropy changes 


across the stator and referred to station 2. The term 


through the rotor. The entropy increase through the rotor is com- 
puted by using the corresponding values of the rotor in Eq. 41; for 
example,k,, would be used instead of H. 


E 
Neglecting streamline curvature and slope, Eq. 9 can be re- 


written 
Ela. ou et ay 2 ees sinc f x aU. cose SIN 62 
eS LE I=. wie ; : . 
cl Xe es eaae 
2 Zn 


ae alles Cos Sz, C cos Lo dhe li C, He ‘cag ds . 
Yo" Was. Yo Warm Xe YS Was. | Ake (65) 

The discharge angle Yen and its derivative ane are found in the 
Same manner as previously described for the stator discharge angles. 
The method for solving Eq. 65 is similar to that used for Eq. 36 
with the exception that iterations are carried out until correspond- 
ing values of Y change by less than 0.005,or un il 13 iterations have 
been completed, The extra iterations are necessary because there 
may be a slower convergence of the values of Y at station 2. After 
satisfying the equation of motion, values at station 2 corresponding 
to those represented by Eqs. 42-47, are computed, It should be 
noted that where absolute velocity terms are used in Eqs. 42-47, the 
corresponding equations for the rotor will employ relative velocities. 
Therefore M. is the Mach number of the flow relative to the rotating 


2 
rotor blade. 
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Overall continuity is checked at the rotor discharge using 
the same procedure utilized for the stator; however, the reference 
flowrate is increased according to Eq. 25 to account for the flow- 
rate between the blade tips and the surrounding shroud, 

After overall continuity has been satisfied at station 2, 
streamline continuity is checked. In addition to the functions 
performed for the check after the stator, there are certain quanti- 


ties that must be adjusted for streamline relocation. These are 


IHe _fdHe ey 
(5 i) (ae) aXe (Xnew ae (66) 














. ; H 
He Novi . Be ci ; = (Xnew 7 Keve) (67) 
Jan edisiic, Se 
( cl e Jee I a Xe oom ae new” Xn.) (68) 
Xnew 
Org | US ed (69) 


A solution exists when the equations of motion,and overall 
and streamline continuity have been satisfied; however, this solu- 
tion has neglected streamline curvature and streamline slope, To 
obtain a solution which accounts for streamline curvature effects, 
the terms in Eqs. 8 and 9 that contain Sy or OR must be in- 
cluded when solving the equations of motion, 

The streamline displacements SY and OR, shown in Fig. 2, 
are computed for each streamline by 


Yo + Yn 
@ 





Sr = a - (70) 


and 
APE Cou tha (71) 


Equation 70 is based on the assumption that the radii of any stream- 
line at stations O and 2 are not greatly different. Therefore the 
cosine of the angle between the line connecting these points and the 
axis of the machine is approximately equal to unity, The length L, 
also shown in Fig. 2, is taken to be half the distance from 0.1 inch 
ahead of the stator to 0.1 inch after the rotor. 
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After a solution is obtained which accounts for streamline 


P 
curvature effects, the resultant pressure ratio — is 


2 
compared with the ie. specified initially. Ideally tou for 
a 


each streamline will be the same. However, computer solutions for 


this quantity may vary slightly from streamline to streamline. 





Therefore a mass-flow-weighted value of fa , called re ) song 
2 & 


is found by 


4 
a . a Ro ee \f te = Me: 
r, oi ; Fe (iri) Foe E 

=| 


R P 
If the specified ro. and $2) aifter by more than 0.0003, the 

Fo Pe Ww 
Mach number of the flow ahead of the stator is properly adjusted 








and another solution is found. 
eeter the iterations for Fe 


2 
cuantities are determined for each streamline,using 


el aes tis = (U,\u, ~ U. Voz) ay 





= AH = Tt, - Ite 
( Anis Ty, - Tacs 


‘he ideal change in enthalpy Ahis is the isentropic enthalpy 





drop from the total inlet pressure Fi, to the static discharge 


pressure P Equation 76 is used to compute the efficiency of a 


2° 


i 


——-— have been completed, additional 


(72) 


(73) 


(74) 


(75) 


(76) 


single stage turbine because the kinetic energy leaving the rotor 
cannot be utilized. The efficiency defined by Er. 76 will be 
referred to as total-static efficiency. 


Theoretical degree of reaction r* and head coefficient ke are 


riven by 
x his-Rés 
Téo ~ [2é5 i} 
and 


kis = Bes a Z2Cp Fea = Tee's gl 
1 (253 Y (78) 
For turbine performance curves it is desirable to obtain the 
mass-flow-weighted values of efficiency ( de head coefficient 
ke Oe? theoretical degree of reaction (r*) horsepower (HP). and 
moment (O° The last two quantities are found from 


Cle on com 


5O (79) 


(Me). = (HP), (S50) 
a= 


ee, 


an 


(80) 
The mass-flow-weighted AH , called (A H) > as well as as de? 


(ke a? and (r*) are computed using enuations similar to Ec. 72. 


Referred values are obtained, following NASA practice. 


For 9 #= 1.4: 


(Phi (81) 


(Me) — (82) 


RPM por = LoL (83) 


: 
" 
8 

T 


Veep —" (84) 
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_ 7 Teo 
Ts7To. 5/84 
(85) 
& = tan - Po 
Psro. 14.7 (86) 


4, MOD I and MOD II Turbines 

The method of analysis as presented was used to determine 
performance curves for the so-called MOD I and MOD II turbines. 
Both turbines are single stage axial-flow machines. Experimental 
tests were conducted on the MOD II turbine by Commons and Messegee 
and are described in Refs. 4 and 6. The test results are plotted 
with appropriate predicted performance curves, 

The so-called MOD I turbine was designed for free-vortex flow 
and has higtly twisted blades. Outer diameter of this turbine is 
9.898 inches. The hub diameters of the stator inlet and rotor dis- 
charge are 6.930 and 5.970 inches, respectively. The stator con- 
tains 13 blades and the rotor 22 blades. Blades of the MOD I tur- 
bine are generally thin. The stator and rotor profiles used to 
predict outlet angles and loss coefficients are shown in Fig. 4. 

The MOD II turbine is approximately the same size as the MOD je 
but its blading is considerably different. The blades of the 
MOD II turbine are thick with blunt leading edges and constant 
profiles over the blade height. Outer diameters of the MOD II 
turbine stator and rotor are 9./01 and 9.836 inches, respectively. 
The stator has a hvb diameter of 6.796 inches, and the hub diameter 
of the rotor is 6.598 inches. There are 19 stator blades and 18 
rotor blades. Stator and rotor blade profiles for the MOD II tur- 
bine are shown in Fig. 5. 

Throughout the remainder of this thesis the MOD I and MOD II 
turbines will be referred to simply as MOD I and MOD II. 

The minimum throat opening "a'' of the blade channels is a very 
critical dimension. Slight variations of this quantity have a con- 
siderable effect on turbine performance, Since this quantity is so 
sensitive, values of '"a'' measured from the actual hardware were used 


for the analysis rather than those obtained from the manufacturing 
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drawings. Then errors due to manufacturing will not be a factor in 


comparison of predicted and experimental results. Figure 6 shows 


the throat openings "a" as a function of radius for the stator and 
rotor blading of both turbines. 

Predicted stator outlet angle X, as a function of radius vy) 
for the MOD I is plotted in Fig. 7. The assumed linear variation of 
N between the hub and the mean radius and between the mean 
and the tip is readily apparent in this figure. Figure 8 shows the 
predicted variations of \., with Mach number M, for the MOD I. 
Figures 9 and 10 are the corresponding plots for the MOD II. 

Relative discharge angles 7 were computed for two radial tip 
clearances k for each turbine; namely, 0.020 and 0.033 inches for 
the MOD I, and 0.015 and 0.033 inches for the MOD II. The predicted 
flow angles *“ are plotted in the same manner as previously 
described for the flow angles oa, . Figures 11 and 12 show [8 


as a function of radius to and as a function of Mach number Mo» 


respectively, for the MOD I, where M, refers to the Mach number 


of the flow relative to the rotor. Figures 13 and 14 show the 
corresponding plots for the MOD II. The predicted effect of radial 


tip clearance on the discharge angles phe can be seen in Figs. 11 - 
14. 


Stator loss coefficients ae were computed for the radial 


locations corresponding to the hub, mean radius, and tip. Figure 


15 shows ~ as a function of radius vy for both the MOD I and 


MOD II. The straight lines in this figure between the values of 
at the hub and mean radius,and between the mean and the tip, 
reflect the assumed linear variation of “_ in these regions. 
Variation of rotor blade inlet angle A. with radius TY is 
plotted for both turbines in Fig. 16. The difference between the 
untwisted and the free-vortex blades is easily noted in this plot. 
Also shown in this figure are curves representing the variation of 


stalling incidence i. with r The change of the MOD I rotor 


1° 
blade profiles with radius is reflected by a considerable variation 


of i. whereas just the opposite is true for the MOD II. 
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Curves for the MOD I showing predicted rotor loss coefficients 
a as a function of incidence ratio ee for the hub, mean radius, 
and tip are shown in Fig. 17. Since the loss coefficients between the 
mean radius and the tip are dependent on tip clearance, there are two 
curves for the tip. One curve holds for the tip clearance of 0.020 
inches; the other is for the larger tip clearance of 0.033 inches. 
For negative incidence ratios between -0.5 and -2.0 the curves for 
the tip are estimations. This was necessary because the computations 
by the method shown in Appendix B, Section 2, gave unrealistically 
low values of Ie as the flow inlet angle Zz, approached -90°. 

The situation is more easily understood when it is noted that the 
blade inlet angle is -33.7° and the stalling incidence is i i at 

the rotor blade tip. Figure 17 shows that the loss coefficients for 
the hub (1, =3.300 in.) are relatively large. The loss coefficients 
at the hub are larger for most incidence ratios than those at the 
tip for a tip clearance of 0.020 inches. The larger value of 

at the hub reflects the higher losses that are associated with an 
impulse type blade. It can be seen in Fig. 4 that the blade shape 
varies from a reaction type profile at the tip to an impulse type 
profile at the hub. 

Loss coefficients for the MOD II rotor are plotted in Fig. 18 
for tip clearances of 0.015 and 0.033 inches, The predicted 
Similarity of the curves for the hub and mean radius is to be ex- 
pected since the blading differs only in solidity. Although the 
blade profile is the same at all radii, the losses due to tip 
clearance result in larger loss coefficients for the tip profile. 

For convenience, the blade properties used for calculating the 
MOD II rotor loss coefficients were those at the hub, mean and tip 
radii of the rotor discharge. Since the annulus area at the rotor 
discharge is larger than the annulus area at the stator discharge 
and since the flow incidence is a significant parameter for the 
rotor loss coefficients, it would have been more appropriate to use 
the blade characteristics at the hub, mean, and tip radii of the 
rotor inlet. However, the error is insignificant because the blade 


profile does not change along the radius, 


3h) 


It may be noted fom Fig. 4 that the minimum radius,for which 
a blade profile is given,is 3.597 inches whereas the radius at 
the hub of the MOD I stator exit is 3.300 inches. The outlet angle 
and loss coefficient for the hub were found by extrapolation, using 
the values computed for the radii of 4.125 and 3.597 inches and 
assuming a linear variation of these quantities with radius. The 
relative outlet angle for the hub of the MOD I rotor was found in 
a similar manner. 

Performance curves for the two turbines were determined for 
the rotor tip clearances mentioned earlier. The total inlet to 
static discharge pressure ratios investigated were 1.30, 1.40, 
1.50, and 1.60, with the exception that pressure ratios of 1.31 and 
1.51 were used for the MOD II with 0.033 inch rotor tip clearance. 
These pressure ratios agree more closely with those experimentally 
investigated by Commons and Messegee. 

The exial distances L used for the determination of the 
curvatures depend on the axial clearance between the stator and 
rotor as well as on the blade geometries. Axial clearances of 0.4 
and 1.0 inches were used for the analyses of the MOD I and MOD II, 
respectively. 

Curves representing the performance of the MOD I are plotted 
in Figs, 19 through 26. Performance values plotted are mass-flow- 
weighted values unless stated otherwise. Figure 19 shows referred 
flowrate as a function of referred RPM. The increase in flowrate 
due to an increase in rotor tip clearance can be seen in this 
figure. Although the flowrate is greater for the larger tip 
clearance,the toique developed is greater at the smaller tip 
clearance. The decrease in torque for the larger tip clearance 
results from the increased losses and decrease in turning angle of 
the flow through the rotor near the tip. The predicted effect of 
tip clearance on torque is shown in Fig. 20 where referred moment 
is plotted versus referred RPM. 

The variation of total-static efficiency with referred RPM for 
the two tip clearances can be seen in Figs. 21 and 22, The referred 
RPM at which maximum efficiency occurs increases when the total 


inlet to static discharge pressure ratio is increased. Blade losses 
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and the kinetic energy of the flow leaving the rotor affect the total- 
Static efficiency. As pressure ratio is increased the absolute veloc- 
ity leaving the stator and the relative velocity leaving the rotor 
increase. Therefore the peripheral speed of the rotor must increase 
to obtain conditions where the absolute velocity leaving the rotor 

is in an axial direction and where the relative flow ahead of the 
rotor has zero incidence. The RPM where the flow has zero incidence 
on the rotor will not necessarily be that at which the absolute 
velocity leaving the rotor is axial. At any RPM,flow incidence and 
absolute discharge angles vary from streamline to streamline,and the 
above statements refer to mass-flow-weighted values. 

It may be noted in Figs. 21 and 22 that the peak total-static 
efficiency decreases somewhat as the pressure ratio increases, At 
higher pressure ratios the ratio of kinetic energy leaving the rotor 
to the work done on the rotor increases. Since the kinetic energy 
leaving the rotor is lost energy for a single stage turbine, the 
total-static efficiency declines. The effects of pressure ratio and 
tip clearance on efficiency can be seen in Fig. 23 where total-static 
efficiency is plotted as a function of the isentropic head coeffictent. 

The variation of referred power with referred RPM can be seen in 
Fig. 24. Peak power does not occur at the same referred RPM at 
which peak efficiency occurs. The peak referred power occurs at the 
referred RPM where the product of total-static efficiency and referred © 
flowrate is greatest. 

Theoretical degree of reaction is plotted as a function of isen- 
tropic head coefficient in Fig. 25. It may be noted that the theo- 
retical degree of reaction increases with increasing pressure ratio and 
decreases with increasing tip clearance for any given isentropic head 
coefficient. The predicted effects are considerably different from 
the results of the radial turbine tests conducted by nie” 

Riley found that theoretical degree of reaction was independent of 


pressure ratio and axial clearance for radial turbines. 


10, i ley, M. W., The Effect of Axial Clearance on the 
Performance of a Dual Discharge Radial Turbine (USNPG Thesis, 
vecember 1966), p. 70. 
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Performance values for each streamline are obtained from the 
computer solution. However, plots utilizing values for each stream- 
line would be difficult to analyze. The deviation of the hub and 
tip values from that of the mass-flow-weighted average may be seen 
in Fic, 26. In this figure,hub, tip, and mass-flow-weighted values 
of theoretical degree of reaction are plotted as functions of 
referred RPM for a total inlet to static discharge pressure ratio of 
1,40. 

Performance curves for the MOD II corresponding to those pre- 
sented for the MOD I are plotted in Figs. 27 through 39, Additional 
plots have been used for the MOD II because of the inclusion of 
experimental results. An axial clearance of 1.0 inches was used for 
the theoretical prediction. Therefore,only experimental data for 
that axial clearance are shown. 

Comments made concerning the performance curves of the MOD I 
are applicable to the MOD II performance curves also. Differences 
in the performance of the two turbines will be discussed later. 

Plots of the variation of referred flowrate with referred RPM 
are shown in Figs. 27 and 28 for two rotor tip clearances. The 
quantitative values as well as the curve shapes agree well with the 
experimental data. The maximum difference between predicted and 
experimental referred flowrates occurs at a pressure ratio of 1.51 
for a tip clearance of 0.033 inches. There are two experimental 
points for this pressure ratio and tip clearance that differ from 
the predicted curve by about 2 per cent. The predicted and experi- 
mental values for all pressure ratios and tip clearances have an 
average difference of less than 1 per cent. 

Figures 29 and 30 show curves of referred moment versus 
referred RPM. The trends expressed by the predicted curves <re in 
excellent agreement with the experimental data. Although the quanti- 
tative agreement between theoretical and experimental values is very 
good for three of the curves, the experimental torque is generally 
lower than the predicted torque. The average difference between pre- 


dicted and experimental values is about 3 per cent. 
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Figures 31 through 34 show total-static efficiency as a func- 
tion of referred RPM. The shapes of the predicted curves generally 
agree well with the experimental results. However, there is an indica- 
tion that experimental efficiencies decrease more rapidly at high RPM 
than is predicted by the theoretical curves. In the high RPM region 
the upper part of the rotor blade has a large negative flow incidence. 


In the prediction analysis when the incidence ratio — had a value 
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less than -2.0, the value of -2.0 was used for computing loss coeffi- 
cients. This limitation may be the reason that the predicted effi- 
ciencies in the high RPM range do not decrease as rapidly as the test 
data indicate. 

The quantitative agreement between predicted and experimental 
efficiencies varies considerably between different pressure ratios 
and tip clearances. There are two data points at a pressure ratio 
of 1.31 and tip clearance of 0.033 inches where the experimental 
efficiencies are over five points below the predicted values. At a 
pressure ratio of 1.50 and a tip clearance of 0.015 inches the 
average difference between predicted and experimental efficiencies 
is 1.5 points. Giving equal weight to all experimental values the 
average difference between experimental and predicted values is 
2.6 points. 

The calculations gave a decrease in efficiency by about two 
points as the tip clearance was increased from 0.015 to 0.033 
inches. The decrease in experimental efficiencies for the increased 
tip clearance varied with the different pressure ratios. However, 
the average decrease in efficiency is close to the predicted decrease. 

Total-static efficiency as a function of isentropic head co- 
efficient Kes is shown in Fig. 35 for pressure ratios of 1.40 and 


1.60. Figure 38 shows degree of reaction r* versus k No experi- 


mental data are plotted in these figures because S Fimental values 
of mass-flow-weighted r* and Kis were not available. 

Plots of referred power as a function of referred RPM are shown 
in Figs. 36 and 37. The comments made earlier concerning the referred 
moment plots apply to these curves also,since the turbine power is 


proportional to the product of torque and RPM. 
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Hub, tip, and mass-flow-weighted values of theoretical degree 
of reaction are plotted in Fig. 39 for a pressure ratio of 1.40. 

The experimental values for the hub and tip er2 also plotted in this 
figure. The trend of the experimental data is in close agreement 
with the predicted trend; however, the quantitative agreement is poor, 
especially at the hub, The decrease in degree of reaction,as the tip 
clearance is increased,is considerably greater than predicted. Both 
predicted and experimental results indicate that an axial flow tur- 
bine diffeis from a radial turbine inasmuch as the theoretical degree 
of reaction changes if rotor tip clearance is changed. 

Figure 40 shows the absolute flow angles “\, and > as a 
function of radius for the stator and rotor outlets, respectively, 
of the MOD II. Experimental data were obtained for the same opera- 
ting conditions; namely, Beggee: toes inches, RPM/{S =13,394. 
The experimental and predicted stator outlet angles have an average 
difference of less than half a degree over the blade height. Pre- 
dicted and experimental agreement for the absolute flow angles at the 
rotor outlet is not as close. If the experimental point near the hub 
is neglected, the average difference between predicted and experimen- 
tal values of 3 is about 3.5 degrees. The experimental flow 
angles X- are larger than the predicted angles over most of the 
blade height, which explains why the predicted torque was generally 
larger than the measured torque, 

Velocity distributions at the stator and rotor exits of the 
MOD II are shown in Fig. 41 for the operating conditions listed in 
the preceding paragraph. The experimentally found stator discharge 
velocities varied only slightly more from hub to tip than predicted. 
The average difference between predicted and experimental stator exit 
velocities was about 1 1/2 per cent. 

The experimentally determined absolute rotor exit velocities, 
which are shown in Fig. 41, are considerably larger than the pre- 
dicted velocities except at the hub and tip. Although experimental 
points could be taken only at distances greater than 0.16 inches from 
the inner diameter of the shroud, the trend of the data points indicates 


that the experimental velocities are less than predicted at the tip. 
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Mear the hub the experimental velocities were less than predicted. 
the decrease in experimental velocities near the hub and tip is 

the result of separated flow in these regions. The reader is 
referred to kef. 6 for additional information concerning the experi- 
mental data and for photographs of the rotor showing indications of 
separation. 

Experimental relative rotor discharge angles Lx were deter- 
mined for the known rotor speed from the measured values of absolute 
rotor outlet angles and velocities, wnich are plotted in Figs. 40 
and 41, respectively, The experimental and predicted values of <<. 
are shown in Fig. 42. The magnitudes of the experimental values 
of .~“. near the hub are larger than predicted. Over approximately 
the outer three fourths of the rotor blade height the magnitudes of 
the experimental values of yes are less than predicted, Ata 
radius of 4.4 inches, for instance, the experimental value of 
is about -57 degrees whereas the predicted value is about -67 
degrees, The average difference between experimental and predicted 
values of < is about 6.5 degrees, The lower values of * pre- 
dicted near the hub,provide some compensation in the overall performance 
for the high values of “a predicted for the outer part of the blade, 

The MOL LI was not designed to achieve the highest possible 
efficiency, The objectives of the design were to investigate the 
effects of b unt untwisted bladings of constant profile. This type 
of blading kas the following advantages: 

1. Blade cooling passages are easily accommodated, 

2. aA wider operating RPM range is posible for a specified tur- 
bine efficiency variation. 

3. Constant profile blades can be manufactured more economi- 
cally than twisted blades, especially if exotic high strength materials 
are needed for elevated temperatures. 

The following discussion concerning comparison of the MOD I and 
MOD II is based entirely on the predicted performance of these tur- 
bines using the methods described in this thesis. The reader should 
be aware that Ainley's methods for predicting loss coefficients were 


not developed for use with rotor blades having blunt leading edges. 
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Therefore the accuracy of the rotor loss coefficients precicted 
for the MOD II is probably not as good as for the MOD I. 

Figure 43 shows the predicted axial velocity ratios as a 
function of radius ratio at the stator exit and at the rotor exit 
for the MOD I and the MOD II. The curves are for the conditions 
that occur at the maximum efficiency of each turbine. The axial 
velocities of the MOD I are nearly constant over the blade height at 
the stator exit and at the rotor exit. The blading of this turbine 
was originally designed for free-vortex flow by assuning uniform 
axial velocity components from hub to tip. From Fig. 43 it can be 
seen that this condition is only approximately satisfied since loss 
variations and curvature influences produce slight variations from 
the assumed distribution. Axial velocities for the MOD II decrease 
from the hub to the tip at the stator exit,and increase from the hub 
to the tip at the rotor exit,since the bladings of this turbine have 
constant profiles. The variation of axial velocities for the MOD II 


is more pronounced at the rotor exit than at the stator exit. At 
WA = 


7 varies from 0.735 at the hub to 1.195 at the 





the rotor exit 
: “AS a 
ip, 

The maximum predicted total-static efficiency of the MOD I at 
a tip clearance of 0.033 inches is about 84 per cent. The MOD II 
has a maximum predicted efficiency at that tip clearance of about 80 
per cent, The difference in predicted peak efficiencies is dve to 
different factors. The stator loss coefficients are highe for the 
MOD I than for the MOD II, whereas the loss coefficients for the 
MOL II rotor at zero incidence are somewhat larger than those for 
the MOD I at zero incidence. An accurate comparison of loss co- 
efficients is more involved for the rotor than for the stator. 
Because of its tw.sted blades, the MOD I rotor has nearly zero inci- 
dence at all points along the blade height in the vicinity of the 
optimum efficiency. Only at one radius of the MOD II rotor blade, 
however, is the incidence zero. The incidence angles at larger 
radii are negative and at smaller radii they are positive. There- 
fore part of the blade is always operating at a loss coefficient 


larger than the minimum. Of probably even greater significance 
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for the efficiency decrease is the kinetic energy of the gas leaving 
the turbine, Minimum kinetic energy is lost when the flow is axial 
in direction. The design of the MOD I is such that at all points 
along the blade height the absolute velocity at the rotor exit is 
nearly axial in the vicinity of the optimum efficiency. The MOD II 
has only a small radial portion of the rotor blade where the dis- 
charge angle Xz is zero, At radii greater than the radius where 
Xo is zero, the absolute flow discharge angle is positive and at 
smaller radii it is negative. Therefore the kinetic energy leaving 
the rotor of the MOD TI is ,reater than that of the MOD I when both 
turbines are operating at peak efficiency. 

A comparison of the total-static efficiency versus isentropic 
head coefficient curves for the MOD I (Fig. 23) and MOD II (Fig. 35), 


shows a larger variation of efficiency with head coefficient for the 





MOD I. For example, at k=0.033 inches and os =1.60, the MOD I 
efficiency decreases about 11.5 points as the head coefficient in- 
creases from 3 to 7 whereas the MOD II efficiency decreases only about 
7 points for the same change in head coefficient. This difference in 
efficiency variation results from the effects of blade twist mentioned 
earlier. 

To investigate the influence of the streamline curvatures on the 
predicted performance of the MOD I and MOD II the presented flow 
equations were solved for an axial length of the bladings L of 9 x 10° 
inches. Increasing L to this value has the same effect as neglecting 
streamline curvature effects as can be seen from Eqs. 124 and 125. 
Differences between performance values neglecting streamline curvature, 


and those where curvature effects were accounted for, were less than 


0.2 per cent for both turbines, 
5. Conclusions and Recommendations 
Effects of streamline curvature were found to be insignificant 


for the MOD I and MOD II. The small effect on predicted turbine per- 
formance due to streamline curvature indicates that the exact value 
assumed for the curvature factor K will not greatly affect predicted 
performance, The results of this analysis also indicate that Vavra's 
method of approximating streamline curvature is of sufficient accuracy 


for methods of analysis where the flow equations are satisfied at 


Stations between blade rows. 
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It is recommended that the MOD I be tested. Results of that 
investigation would provide additional information concerning the 
general applicability and accuracy of the method of analysis pro- 
posed in this thesis. 

Predicted and experimental flowrates for the MOD II were in 
close agreement. Since the restriction factor g is a significant 
factor in the predicted flowrates, the experimental results verify 
the validity of the theory used in the developement of & . 

Experimental results for the MOD II showed that the predicted 
rotor torque and turbine efficiency were generally 2 to 3 per cent 
too high. The angles and velocities predicted for the flow at the 
Stator outlet were in excellent agreement with the experimental re- 
sults. Therefore the high values predicted for rotor torque and tur- 
bine efficiency must be tied to the rotor solution. The measured 
values of outlet angles and velocities for the rotor discharge also 
indicate that the predicted rotor solution is not completely correct. 

Exactly what parameters are in error in the predicted rotor 
solution, and to what extent, cannot be stated without more experi- 
mental data. More traverses should be taken at the stator and rotor 
outlets using recently calibrated probes. With this information,loss 
coefficients and relative discharge angles for the rotor could be 
determined on a streamline basis. Based on the experimental results 
included in this thesis,it is suspected that the magnitudes of the 
predicted discharge angles Be are too large. If the predicted 
turning angles of the flow through the rotor were less, the predicted 
torque and efficiency would decrease, Also, the separated flow at 
the hub and tip at the discharge of the MOD II rotor indicates that 
predicted loss coefficients should be higher for the streamlines at 
these locations. 

To accurately predict turbine performance at off-design as well 
as design conditions,a streamline analysis is necessary. The experi- 
mental verification obtained for the proposed method indicates that 


the method possesses much potential and that additional developement 
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is warranted. For future developement work on this type of analysis 
it is recommended that the following changes be considered: 

1. Use seven streamlines instead of five. 

2. Calculate discharge angles - using the methods of this 
thesis but neglecting blade curvature effects. That would decrease 
the angles by 2 to 3 degrees and be in closer agreement with experi- 
mental results. 

3. Introduce a multiplying factor for the rotor loss coeffici- 
ents which would more accurately account for the separated flow at 
the hub and tip but would not greatly change the average rotor loss 
coefficient. For example, rotor loss coefficients would first be 
calculated using the methods of this thesis,and then the loss co- 
efficients for streamlines 1, 2, 6, and 7 would be multiplied by, 
say, 1.4 and those for streamlines 3, 4, and 5 would be multiplied 
by 0.6 or 0.7. 

4. Calculate rotor loss coefficients for an range of 
-3.0 to 2.0 instead of -2.0 to 1.6. The increased --~- range would 
improve performance prediction for turbines with untwisted rotor: 


blades at off-design conditions. 
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APPENDIX A 


DEVELOPEMENT OF EQUATIONS 


1. Equation of Motions 


The equation of motion for relative flow is 


VRR=eWxly x +2H)4+ TVS 


In cylindrical coordinates: 














HY = 2 Oy + ¥i 3 H ic = + or 
a if QF = a Cc Sy 
ay! & aN yo 
i — j : A 
— | i 2) vy ale El 
my © 3S 22 | aaa 99 {h 
3 jo) ay 1 | r1oWe, trv vie 
a C2 \\ cto - a) - Wo = ean =~ oe al 
~ Vise C | ech) s2 1) 
; ; \ : : = 
oe ho - a(eWs) _ 9Wr | _ yw xe awe]? 
a 0 ary ss a 96 Aiy2 or J] 
WV x 2 DW co 2 Ws) oe ig (2w Wr) 
> i ee 7 as “ 9S ie ie : 


Equating Eq. 8/7 to the components in Eqs. 88-90, 
_ } 7 ' Nea Ai hy ; aS 
a a atte. Yin 2M adevs)] We | DOV) WV yy, 12 
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Eckert, op. cit., pp. 149-155. 
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(88) 


(89) 


(90) 


(91) 


(92) 


(93) 


















(94) 


For assumed axisymmetric flow, ot = 0, 
We alr.) _Wr rWe) 
lg. O= a v or -2WW, 
mS 2 R_ wr oVii6 a Ws arWe) ~ 9s 95 
at Se We “Se Woiemenc as PSS ey 
+ SHe. Wo otrW Wie _ CEE 
i owe te cou) -Wa SEP + Wp 8 + 2 Wo + T Se (96) 
Equation 94 can be written 
Siew.) Wr ote W Eee 
a We a -2wYr Wa (97) 
Replacing an” in Eq. 95 by its equivalent in Eq. 97 
Pa Le 
3 ate = W, a“ -\' Wn a Wa, Gan Wr aus v D2 toes VV. gt (98) 
oF ——!|COUY VVA of 
or 
oHe — Wolsc adv W)) hy aule = Vee gt (99) 
oF en a), oT +VWVr “Wr ys JV gs 
Multiplying Eq. 96 by We and = 99 by W me 
a 2 Nee Z Wi We alr val Wi. Nr >We =p Wr 2M wl. Wi, +WeT = | (100) 
Sir ¥ af 
Wa Se = Wie 2/2) ays, yy, Me oe Wir s¥igT 3 (101) 
Adding Eqs. 100 and 101 
st Vip EE = (We BS ~ Wa =) (102) 
en a 
For adiabatic flow Hp is constant along a streamline, therefore, 
(103) 
(104) 


Vi coy Ne = 


\\\ aa -V 4. se 


i 





For axisymmetric flow, 
= io SN 


>t | 
(105) 
(106) 


From Eq. 104 
ate 2) We ote 
ee Vip =e 
Substituting Eq. 105 into Eq. 102 
, dHe Wr He \ ; 9 3S° 
Ne Se “We (- in eo Ss: + We Sa 
or 
os a We Pi (107) 
enone 
90 


= 
2 


Using Eqs. 105 and 107, Eq. 99 becomes 


r_ SHR Vive otev' } ~Wr Nas -\Vr ave VA -2u Wolly ve Me T= (ee) 

















Wao ss CVVg ere a2 of Via ‘fh Wa A oY 
Multiplying Eq. 108 by { “e \ , gives 
SHe _ We oth) yy, AWE yy Ove 2 pW, eT 2 109 
ae Sa ao” oc es ae 


Eq. 109 is identical with Eq. 96,therefore Eq. 109 is the equation 





[ ot oes 4 (We. e) ‘ 
that must be solved. With | Mp =~ — -|,Eq. 109 becomes 
a se om, “> \/ a Wien aie -93S 
eel. - Ni - ae 2hW gow, -z OT S20 (110) 


of lie: ae eal 


Rewriting Eq. 110 for the rotor discharge and using the substitutions 


{ \2 = He Ws 
4 =H & ‘eS. rag i a 
KR. H. 3) = d wee Ce ace 





We Ne Se, wee ue Wee 
a —o bor 3 Wee ae Wa) ~ 4dwWy, -2 5 2. eS +2{-2 -—-|— 9 

















ies ed al gz Ee a ae cP 2p a, (111) 
The last term of Eq. 111 can be expressed by 
te (yj | ail tains Wea | (112) 
With a ’ ) ape 112 becomes 
Ve me oe ' 
' | : Was 
ee His =\ir oo =< = \2 es a ee Az. 
Cpt re We Cel™ He Wo. Cp Gee tie (113) 
Using Eq. 113, and after rearranging, Eq. 111 becomes 
S(Vib } ov Was aSe2 a WV Vite at 2 W5, 
Sie ee 32 St 
3 ; ~F ~Y> fle ] 2 3 {> 
a ee “35 
+4. — ae Seen 114 
ta: eS Bs ele He SNS Sv © ( ) 
Noting that 
y GUE ; a! : a) : Lae 
[oa oe ~ oe (115) 
Eq. 114 is : 
a fey 
iV ie yy Lf V4 Jae Waez (Wu. ie Wu- 
Si - 37 Cp Cos Ao aT, = (2 Qt 
Ome ee oz 
-~ 72. ~ a7 if a - ZA \ J ai = O (116) 


BE 


a Vise 
This equation is made non-dimensional by multiplying bY. 
[Byn 
where the subscript m refers to the mean streamline. In the sub- 


sequent derivations the subscript 2 will be omitted. Then 











‘2\ x ie os er aN seer 
Yow >(We te Wa r dvr _ Wr, View cs nee Vivu es CN. Wan 
Win 2 We, OB Wal, Cp cesta oY 9 Wan fy CL 
ee =e See | 2H Wo 3S 
+ AS ee aman — 2 ts = a ee a eee eo ag ——- ee a (117) 
Was. 1 Wat. Vien =e { ie 
Now let WW 
— -- vA 
Y= Wer, (118) 
2 aus 
\ Sara (119) 
ate ae (120) 
cp 








CY) gp Voy Me YE ast yr? aXytané) 
: oe id x aX% 





YY Tene ce ie i ae se 
ye 2 Sie pete yeran [2220 (121) 
WA. Wi. Vb nn aK 


The fourth term in Eq. 121 can be expanded and rearranged as 


2 ¥ TE | 3 -(KY Tene) = za ee vba A tany| 


yc 
_3TANS ail ap 
=? eee Pi gis AN +—— TAN (122) 
y 23 WK a ee 2 Neal yg x F 
with 
oP ks IF 
>» me 2% 
and 
2 2) ,Ge 
SN Se a = TAN £ é 
aN | = ox TAN , SX 
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~ gee 122 then is 


dWr_ _Y* ds* 2TwNG VB. 
a, ranté)-2, Se ace 5% aay ae s 2? safe 


seat a ave , are |S -Y rp mies 0 (122a) 


Multiplying by (- i ) and with |+TAN “b= ~ 


ae a(Y*) MEY Sa as"). Tanp 22 + a sin°f6 
¢ A J 


Av 








tn sin Bees | 2 Um cos B 2c He " 2 He cosiS ey zg)9S" (123) 
PSE To Vee CR ce dX 


The terms cos‘) and i Mr) represent the effects due to 
streamline curvature and can be approximated by 
_j SWroe , Sh] Kes (124) 
= = : 
Wp 32 is 
2 
1 : 
2 = ne ee 
CO = ae, (125) 
L*+($*) 


If the streamline curvature in meridional planes is zero, the terms 


represented by Eqs. 124 and 125 will take on values of zero and one, 
respectively. With Sr positive as shown in Fig. 2,the plus sign 
is used for K at station 2 and the minus sign at station 1. The 
streamline slope as represented by Eq. 125 is the same at stations 
1 and 2 since the streamline pattern is. assumed to repeat itself after 
station 2. 

Using Eqs. 124 and 125, Eq. 114 then becomes 


1. tlle) aml ow 


AUS! NBcosk 2 Un Ucos?B 2Cos 2 Contam de Peasant He cost Benn a 


Wan Y Wa ¥ ” Wat V2 aX | Wiy* G2) 
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For the enthalpy terms to be in terms of BTU/1b_, they must be 
divided by 2gJ to keep the equation dimensionless. With 
C, = 2g J, Eq. 126 is 


a % 2 § LEY ds*|_ dé a is 
ee 2 -cosie enn )-( - jas 2TANB ay SINE 


- 4 Um SINS 6 cos B oe ZU me cos B | & cos - 5 dHe ’ He cos 4 eX 
Way. Y Wan ¥~ WamY? dX 


Noting the corresponding terms for the absolute flow of the 
stator, and with U=0, the equation for station 1 is 


a as = costal Karat) {So | 2 -Z2 TANS a8 





S Ce oe 
2 2 C, Cos So dH C, Heosta c ds* 
-—SIN® ¢ ———— + Oe sin (128) 
X evn. dX aes 


2. Flow Function } 


Flowrate in 1b _/sec for the expansion process shown in 
Fig. 44 can be written 


we PAV 


Fa a ~\ 


—_—_—-> W 


(129) 


(PLENUM) 


\ 
‘ 
L} 
t 


A,P, We 


Expansion From Plenum 
Fig. 44 
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Assuming an isentropic expansion, 


~ ta 
PT ee \ ; 
| P| 


and, 


; ‘ x-\ 
es FU! ae 
ae L a ai 


then using the thermodynamic relations 


ae 
aa oe ae ana Cp hata 


The density ie is found with the isentropic assumption by 


ie ‘c 


or 


Using Eqs. 133 and 135, Eq. 129 becomes 


~ 


ae 
We 


\ 
This can be rearranged to give 


or 
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(EY Jeane whi-lé 


(ged FO B\S_ _P\ 
N = ie ay he KR me 


(130) 


(131) 


(132) 


(133) 


(134) 


(135) 


(136) 


(137) 


(138) 


3. Restriction Factor Ee We 
The expansion process through a blade row isa shown in Fig. 45 
(a). The letter u represents the velocity at some point within the 


boundary layer and ve represents the velocity for an isentropic 


h : 
expansion. The various dimensions, angles, and velocities for the 
throat location are shown in Fig. 45(b). 


af 


(7) 





(a) 


Conditions at Exit of Blade Row 
Fig. 45 


The mass flowrate per unit blade height can be written 
° V t Z% : d (139) j a 
i Pan Uy costal Sahay - ER] +E fury 

oO 


Ve = density at P, 
a = density corresponding to the velocity u (ilGas/ft’) 


where: 


for isentropic expansion (slugs/£t>) 


X4y = blade angle 
$ = boundary layer thickness 


12 avra, M. H., Problems of Fluid Mechanics in Radial Turbo- 
machines (Rhode-Saint-Genese, Belgium: Von Karman Institutéecfor 
Fluid Dynamics, 1965) VKI Course Note 55b, pp. G46-50, 
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For Mach numbers of 0.8 or less, the discharge angle can be 


calculated quite accurately by 
ee 
CosKy 


a = minimum throat width 


t = blade thickness 


then with / = — and using Eq. 140 in Eq. 139, 


u 


hs ) Geos a S Cos M4 - -t Scos oc t a 
ge" (ta ens Scocad a Cos KT a ape =) Z ferMen® Go acl 


which can be expressed as 


; | $ o 
i U, a Ea (-fZ, 71) 


Assuming constant pressure at the throat, 


Be 
{PT 
and 
fo, lt 
fer ov 
where ye Ven 2 Y 
ir. to 29J¢p to 2g3¢ ee = 





In Eq. 144 the term Xx, 


Ten oe 
Ke=l- = Pe 


The boundary layer is assumed to be turbulent and the profile can 


therefore be expressed by 


Vee We 
in CE) =F 
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Too Te Vth ee 
te agice ‘ Tee ee a |-Ke ( 
s 


u \e 
i 


f 





(140) 


(141) 


(142) 


(143) 


(144) 


(145) 


(146) 


Then with the displacement thickness of the boundary layer defined as 


ses fie (oa/ toma ‘ 


Eq. 142 can be written 
¥* 





é Zo 
AAS: Veh a|1- a (148) 
The loss coefficient to the throat is 
2 
mel = Na (149) 
ies 


The loss coefficient can also be expressed in terms of the kinetic 


energy lost, by 


ae = 
E my 150 
s Ve 2 | . ve ( ) 


: ae 
where E represents energy rate at the discharge due to the average 


velocity Va or 


Eh Ve la- 45] Mek 2 £fupt ") dy (151) 


Equation 151 can be rewritten to give 


3 1 3m 
: Vth $ Y 2 
= eo | \- = | = l-Xe (152) 
E fen a z Q ( ) (\-Xe 7") ( 
Oo 
Then using the energy thickness which is defined as 


sme sli- (ufo, an| (153) 


the loss coefficient can be expressed by 
v. oe v2 ee 

zi A, ee ali-2 a ] te a —ta £2 
m —sa- Van =| 
= 2 Vere af 2 


of = le lo Se (154) 
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The restriction factor — represents that part of the throat 
opening in which would occur the uniform theoretical velocity, 


therefore 


on 
— iti 1 
Ea Bie (155) 
Now defining an energy parameter H*** by 


KKK qs 
H = = (156) 


Eq. 155 becomes 
pesm oa 


= Lens ee ars 
The loss coefficient “4% in Eq. 157 accounts for the losses that 





(157) 


occur from the inlet to the throat of the blade channel. No means 
exists for predicting this loss coefficient. Half the total loss 
coefficient for the blade row tt will be used to represent these 


losses. Eq. 157 then becomes 


*** 
E= oad (158) 
H**E \ + Se 
mene 
4. Method of Evaluating H*** 
By use of the binomial theorem, the denominator of the integral 


part of Ecs, 147 and 153 can be expanded,yielding 


Bm 


=| 
Zee am 2 4m 3 Gm 4 
(Ke 7) dake ake ake 1 Ke pe = 
The integral of Eq. 153 is then 


Nt? keys ea " Keg +Xe" ) ae | dy 


Then, integrating and evaluating gives 




















om 2 3 4 
eS ca ower, Ken ee Ker ek 
1-Xeg m+]  Sm+l  Pmtl Dmtl ime 
O 
Equation 153 can now be expressed as 
kk % 3 4 3 
pee ef eter he ee ke Xe. Kee 
2 i. Sm+l Smell Fmel Dmel \lme) Zm+] Smt) Tme!] 


Ie Meteet , @p. cit., pp. to9=160; 
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x * 
2 3 4 
3 =\+(Xe-N)[+ re, KeF _ 4 pee ecipeanes 
3 Sms} Smt] Jmsel Det! Ime 
also, 
x 2 7 4 
3" a X -|) a Ke + Ke . Ke + Ke does 
. © Re-l Bm +l Smal 7+) Qe) me! 
In a similar manner, the expression for Eq. 147 is 
* 2 2 4 
See oly = Xe | Ne) Nema 
$ Rs Xe-| mi Bme+l Smell Ame! Sm-+) 
The energy parameter can now be expressed by 
ye eee ae 
eee Xo.) (Guage Sa ee a oe 
al = e-! sme] Sm+l Amel Qmr) mn + | 
_ ae 159 
ee et ee et 2 
Kea m+| 3 m+ Sm) Tm Sm) 


The number of terms in the numerator and denominator of Eq. 159 

is considered sufficient to give good convergence for H***, 
Although the exponent m for the turbulent boundary layer is 

dependent on Reynolds number, it is taken to be a constant for 


this analysis. 
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APPENDIX B 
COMPUTATION OF OUTLET ANGLES AND LOSS COEFFICIENTS 


1. Outlet Angles 

The absolute discharge angles for the stator and the relative 
discharge angles for the rotor are computed by using the same 
methods. For these calculations, the assumption is made that the 
outlet angles are not influenced by the flow incidence angles. 

The factors which do affect the discharge angles are: 

1. Blade geometry; this results in the angles being a func- 
tion of the radius since spacing and possibly profiles change from 
hub to tip. 

2. Radial tip clearance; the effects of tip clearance are 
assumed to influence the rotor flow from the mean streamline out 
to the tip,with the largest effect being near the tip. 

3. Exit Mach number; in accordance with the experimental 
results surveyed by Ainley, the outlet angles are a function of 
Mach number. Values of the angles are calculated for Mach numbers 
M of 0.5 and 1.0. Then a smooth curve of outlet angle versus M 
is drawn between these points with an inflection point at M=0.75. 
Below M=0.5, the flow angles are assumed to be equal to the value 
computed for M=0,5. 

Vavra's formula is used for the first approximation of the 
stator outlet angles for M.50.5. 


a *= cos'(FE) (160) 
Ke 


3.3 
sae = ee Q 
K, =Pl\ = io ( 120) S (161) 


The effects of blade curvature are next taken into account 


by using the method given by Ainley, \4- 
= Ay A ew (ue?) 


14, snley, op. cit., pp. 3-4, 
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In Eq. 162, e is the mean radiu: of curvature of the upper surface 
of the blade profile between the passage throat and the trailing 


edge, This quantity is approximated by 
2 
J 
4 = 163 
ab (163) 


where j and b are shown in Fig. 46. 





Blade Geometry 
Fig. 46 


The discharge angle for M=1.0 is found in a manner similar 


to that described by Ainley, 


-1f/ a 
=a > 
Han, = 1.0) ou & =a ) 
The change in & due to the increase in local Mach number from 
0.5 to 1.0 is assumed to be constant along the blade height and 


equal to that computed for the mean radius, 


§ Hp = SHH = Xr Geo.ey ~ me (etd (165) 


Relative discharge angles for the hub and mean radii of the 
rotor are found by the same method as used for the stator. The 
effects of tip clearance are accounted for in the computation of 
ee for the tip position. With Ainley's formula, the flow 
angle yo for the tip is 


f2=TA N- It -X (EEE: rs b. * tt eae raed (166) 
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The new terms in Eq. 166 are: 


h - blade height 
UI 
427 discharge angle before ti; clearance is accounted for 


X - factor which depends on the type of shroud 


Ainley suggests that X be set equal to 1.35 if no shroud is 
arranged at the rotor tips. However, Ainley's value has been 
determined with a mean streamline analysis. To localize the 
effects of tip clearance, X has been changed to 5.4, a value which 
gives about the same overall tip clearance effect as if the value 
of X=1.35 were used for the mean streamline. Fig. 47 shows that 
the areas under the two lines representing X as a function of radius 
are approximately equal. The slope of the solid line representing 
X for this method is obtained because of the assumed linear varia- 


tion of outlet angle with radius between the mean radius and the tip. 


USED FOR MEA 
5.4 S.L. ANALYSIS 


1.35 





Yrup Ton Yuu 


Tip Clearance Factor X as 
Function of Radius 


Fig. 47 


Sample calculations for stator outlet angles are shown in 
Table I. A sample of rotor outlet angle computations can be seen 
in Table II. 
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TABLE I 
SAMPLE CALCULATIONS FOR STATOR OUTLET ANGLES 


MOD I turbine 











M,= 0.5 

radius - in. 3.597 G25 4.950 

spacing (s) - in. 1.738 1.994 Zane 

throat opening (a) - in. 0.455 0.590 0.804 

trailing edge thickness (t) - in, -------- 0.044 ------- 

= | 0.2619 0.2960 0.3360 

3 = 

(E100 A 21.4 13.6 7.48 

K,=1-(0.0027) (A) ©) 0.9849 0.9891 0.9932 

— 

Ky 0.2658 0.2991 0.3384 

rn coe 7an6 TONG ee One 

SKy 

4 () 0.703 0.9106 1.37 

a mar 42) ise” ~~ ees em 
M,=1.0 


Use 4,125 in. from centerline as representative radius. 
Projected trailing edge thickness (t*)-in. 0.13 

a 6 
rearey, - deg. TAGS 
2A 73.5 - 71.5 = 250 





rr ~ cos”! ( 


oO Bi Oem eS 69.6 
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VABLE ii 
SAMPLE CALCULATIONS FOR ROTOR OUTLET ANGLES 


MOD II turbine tip clearance =0,033 in. 


M £ 0.5 

radius — in. 3.3235 4.197 4.918 
spacing (s) - in. 1.595 1.4643 1.7158 
throat opening (a) - in. 0.382 0.560 0.721 
trailing edge thickness (t) - in, ------- 0.048 ------- 
5 0.3295 0.3824 0.4202 
(E) 00) 4pld97 3.278 arene 
J 109.0 50.3 29.8 
K.= |- (0.0027) (a)(Z) 0.9030 0.9481 0.9662 
ar¢ 0.3649 0.4033 oe 

m5 <7! 4) Bepes 1 e662 = 64,2 
A (-2) Wars <0 2.37 3.04 
[3:8 "- 4(2) . 270.3° 368.6 -67.2° 
k - in. 0.033 
k /h 0.0207 
X=(4) (1.35)*type shroud factor 5.4 
cos f,, $.=30° 0.86603 
cos, , 4, =-67.2 0.38268 
(x)(4) aS =6 0.25297 
tan i =-2.4142 
(1-B)tan -1.8035 
tan So 0.57735 
(1-B)tanZ°+B tanf, =C -1.6574 
/o=tan*(c) __-58.9° 
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TABLE II (continued) 
SAMPLE CALCULATIONS FOR ROTOR ANGLE 
Mj=1.0 
Use radius = 4,197 in. as representative radius. 
Projected trailing edge thickness (t*)-in. 
7 -l\/_a 
(ae a) 


06 =-68.6+65.6= -3.0 





radius - in. 373235 “421197 


bz -67,3 -65,6° 
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0,108 
-65.6° 


4.918 
-62.6- 


2. Loss Coefficients 

Loss coefficients were predicted using Ainley's methods. 
The techniques used are completely described in Refs. 1 and 2. 
The assumptions, definitions, and equations which are necessary 
for a basic understanding of the method are described below, 

It is assumed that loss coefficients are independent of Mach 
number and, for a given glade geometry, are a function of the flow 
incidence on the blade. Losses are divided into profile losses, 
secondary losses, and losses due to tip clearance. Mixing losses 
are not accounted for separately but are included in profile and 
secondary losses. 

A parameter Y is defined by 

Y = Loss in total pressure 

Total pressure at outlet - outlet static pressure (167) 
For incompressible flow, which is assumed for these calculations, 


the loss coefficient is 


a = Suen (168) 
1+ ZY; 

The subscript ; on Y in Eq. 168 indicates the type loss repre- 

sented; i.e., a would be the profile loss. 


For zero incidence, the Y for profile losses is 


; Ps 
Yves” Yoaney‘(a) aay reroll $4) - a 


where: 





A, - blade discharge angle not accounting for tip clearance 
th - maximum blade thickness 
c - blade chord 
Equation 169 is applicable to both the stator and rotor. However, 
for stator use the blade inlet angle Ke is zero. The quantities 


e(4,=-fi') 


quantities are plotted as a function of blade solidity. The maximum 


and Ye(gz0) ate taken from Fig. 4 in Ref. 2 where these 


value of ive is set equal to 0.25 for thick blades, 


To determine profile losses for incidences other than zero, Fig. 
6 in Ref, 2 shows Ye versus -— where i is defined as that in- 
Ye (izo) ts. s 


cidence where the losses are twice che minimum losses. 
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Secondary and tip clearance losses are computed by 


bP 6.) * cos*& 
Ye +¥q = [n+ B lea Sas: (170) 


where; 


fm the mean of the gas inlet and exit angles of the blade row 
B - factor representing type shroud, similar to X 


st = 2(tan£,-TANA, )cosBu = Fé) 


(171) 

N = £\ (Ad/p)*//\, 12 

f ( /n:) {0+ 

in Eq. 171: 
é 

A, = (annulus area at discharge) (cos 4. ) 
A, = (annulus area at inlet) (cos fa ) 
I.D, and 0.D, - inner and outer diameters, respectively, at the 


blade exit 


Values of /\ are obtained from Fig. 8 in Ref, 2. Ainley suggests 
a value of 0.5 for B for the type shroud used. With the same rea- 
soning as previously used for the effect of tip clearance on dis- 
charge angles, B is set equal to 2.0 for the losses at the tip and 
zero for the mean radius and hub locations, 

Sample calculations for the prediction of stator loss coeffi- 
cients are shown in Table III. Sample calculations for prediction 
of stalling incidence and rotor loss coefficients are shown in 
Table IV. 
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TABLE IIT 
SAMPLE CALCULATIONS FOR STATOR LOSS COEFFICIENT 


MOD I turbine inlet angle ( &, ) = 0 
radius = 3.597 in. outlet angle ( X, ) = 5Ho 
pitch (s) = 1.738 in. chord (c) = 2.691 in. 


All figures referred to are in Ref. 2. 


ANG e "7, 
Ye Gee) : Pewee i Ce) Yew -o\ Youig ca 


A,=0 Wings (ete 


=0.64S9 
CHORD 2. 2.69) 


Y= [> +8 FH) [SE] ms 


/ (a), Ki 7 #2) | = ]= $f Ges 5.3), , as ro y¢ an a 
A= ¥ (0.0386) d= 0.0086 (Fig, 8) 


Kya = TAN" TARZES) = TAN ({.906) = 62.3° 





& =2| TAN Ke-TANA,|cos%m = 2 [- (TAN -75,3°)|cos 62.3° 
(ae 
= 3.5438 
ai 
Y,= [0.0086][ 3.5436 "estes 7$.3° |- 5,045| 
\ I ] Cos® 62.3° 


vies 16 =Cveo2o 


CORRECTION FACTOR FoR TRAILING EDGE THICKNESS , F=|.03 


(Fig,2) 
(1.03)(0.0929) = 0.0949 


ie 2.6989 


= 0 
\+0.0949 Tee 
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TABLE IV 


SAMPLE CALCULATIONS FOR STALLING INCIDENCE AND 
ROTOR LOSS COEFFICIENTS 


MOD II turbine s = 1.7158 in. So = 30° 
radius = 4,918 k = 0,033 in. fog = -67.2° 
t = 0,728 in. ¢ = 1.967 in. 

max 


blade height(h) = 1.595 in. 
Figures referred to are in Ref. 2. 


BN= c Pofg 
Yoces” ean (B) a wees} t * a 


t _ 0.728 Ao al 
ee ee (CE ee - - = 0,446 =20,1995 
pag, cae ‘ Ox 


= Az 0.154 Fig. + 
; Yeyg-o)* 00365 I62&) (Fig. 4) 








0.446 
Fino) = 30-0865 * (011983)(0.154- 0. o8e5)€ (I, 25) Pe Clee 





> ee 
ae = =) S Fé ; =-70.4 = Di4ZG 
(7 (3 = 0,75) . + Eis ), ROA ; 


Scgvonsy pige-7° (Fig.7) , (,#32%S%7°= 25.5 


Y.+\= |? rele )\| ad [4] B=2.0 oS nell 


cos, = 0.38782 cost*g’/= 0.15017 


N= te Deyo rb. T,D, OUTLET= 3.8235 CN, OD, OoTLET = 4,518 iN, 
A, (+ O.D. ; 


L.D, INLET = 5.3245 in, OD, INLET = 4.8506 tn, 


e 


fo = 7 
Az _ 1(4.918%- 3.32357) cos 70,8° 
A, 


tt (4.8508*-3.3245") cog 30° 


Ne £ (0,148) . ) = 0.0072 (Fig. 8) 


2 
Ne Vy = lo.0072 i, rors 


Con 4. 
+ Vx = (0.008206) |¢ Hen) 





nee. 
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TABLE IV (continued) 


SAMPLE CALCULATIONS FOR ROTOR LOSS COEFFICIENT 


Vi, 

[= (Lis) bs 
Ves = t+ Bo 
TANG, 
TAN: +TAN GE = A 


& 


B= TAN (A) 
cos B,, 


TAN 6,- TANK, = B 
“F7.7 2 (B) cos Bm 
a 
E+/s7] 
Viost fh, 
eae 
eee =(0.003900)| <r ae 
eas Fig. ¢) 
Ye = | Ye/ Yeu=0y] YpG =e) 
i, + Ye + Vx 


=125 


a 


-38.25 
-8.25 


-0.1450 
-1.2619 
-51.6 
0.6212 
75 oS ) 
2012 
7.7016 
4.1726 
0.2345 
oor 
0.1390 


0.3735 


(correction Factor\(ZY)=C° (Fie2) 0.3922 


a 


—_— ss 


rae 


_0,2817 
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0. 


Oo [© 


30 
0.5774 
-0.9008 
-42.0 
0.7430 
2.9563 
4.3928 
19.2964 
2.4385 
0.3434 
1.0 
0.0662 
0.4096 
0.4301 


0.3007 


las 
38.25 
68.25 


2.5065 
0.0638 
3.6 
0.9980 
4.8854 
9.7510 
95.0813 
1.0061 
0.6980 
4.5 
0.2979 
0.9959 


1.0457 





APPENDIX C 


COMPUTER PROGRAM 


A computer program using Fortran IV was written for use with 
the performance analysis. The program was used to predict perfor- 
mance values for the MOD I and MOD II turbines. However, an 
attempt was made to keep the program general. If the methods of 
this thesis are utilized for expressing turbine characteristics, 
this program can be used for other single stage axial turbines, 

Input for the program consists of information representing 
the turbine characteristics, an indicator specifying the detail 
of output desired, and the conditions for which performance values 
are to be obtained. 

Twenty input cards are used to introduce the turbine 
characteristics. The information contained on each of these cards, 
with required dimensions, is listed below. 

1. Stator mean streamline gas outlet angles (ALFAM-radians) 
for mean streamline exit Mach numbers of 0.5, 0.7, 0.75, 0.8 and 
1.0: 

2. Relative rotor mean streamline gas outlet angles (BETAM- 
radians) for the same Mach numbers listed in (1.). 

3. Radii ahead of the stator (RC-in.) for the five stream- 
lines ,such that the flow area is divided into four equal parts. 

4. Assumed radii of the five streamlines at the stator exit 
(RS-in.). 

5. Same as (4.), only for the rotor exit plane (RR-in.). 

6. The predicted stator loss coefficients (ZETAS) for the 
hub, mean radius, and tip. 

7. Ten values of incidence ratio (RINC) ranging from -2.0 
to 1.6,in increments of 0.4, 

8. Ten values of rotor loss coefficients (ZETAR1) for the 
hub corresponding to the incidence ratios of (7.). 

9. Same as (8.), only for the mean radius (ZETAR3). 

10. Same as (8.), only for the tip (ZETARS). 
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11. Length L shown in Fig. 2 (CL-in.) and the curvature 
factor K (CK). 

12. Differences in stator gas outlet angles (DALF-radians) 
from the angle for the mean streamline, for streamlines 1, 3, and 5. 

13. Same as (12.), only for the relative rotor outlet angles 
(DBET). 

14, Ten radii for the stator exit plane, equally spaced or 
approximately equally spaced, ranging from the radius of the hub 
to the radius of the tip (Rl-in.). 

15. Ten values of stator throat opening (Al-in.) corresponding 
to the radii of (14.). 

16. Same as (14.), only for the rotor exit plane (R2-in.). 

17. Ten values of rotor throat opening (A2-in.) corresponding 
to the radii of (16.). 

18, Inlet blade angles for the rotor (BETO-degrees) for the 
hub, mean radius, and tip. 

19, Ten stall incidence angles (STALI-degrees) corresponding 
to the radii in (14.). 

20. The radial tip clearance of the rotor (TIPC-in.), the number 
of stator blades (ZS), and the number of rotor blades (ZR). 

Input card number 21 specifies the output to be printed, and its 
use will be described later, The remaining input specifies the con- 
ditions for which performance values are to be found,and enters the 
estimations of the flow Mach numbers used for the first approxima- 
tions in the iteration process, Input card number 22 specifies the 
number of sets (NSETS) of operating conditions for which solutions 
are to be found. A card containing the following information is 
used for each point specified by NSETS; 

a, Estimated Mach number of the flow ahead of the stator (AMC). 

b. Estimated Mach number of the mean streamline flow after the 
stator (AMS). 

c. Total inlet pressure (PTO-psi.). 

d. Total inlet temperature (TTO-°R). 
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e. Estimated relative Mach number of the mean streamline 
flow at the rotor exit (AMR). 

f. The rotor speed (RPM). 

g. The ratio of total inlet to static discharge pressure (PR). 

There are eight subroutines in addition to the main or execu- 
tive part of the program. The subroutines and their main functions 
are listed below: 

1. Subroutine PARAB is used to determine the coefficients of 
the parabolic equations used to approximate curves. 

2. Subroutine LSQPOL determines by the method of least squares 
the coefficients of the fourth order polynomials used to approxi- 
mate curves. 

3. Subroutine CHAN computes the flowrate and the reference 
flowrate from conditions ahead of the stator. 

4, Subroutine STATOR determines the axial velocity ratios that 
satisfy the equation of motion at station l. 

5. Subroutine FLOWREF computes the reference flowrate at the 
Stator exit and the rotor exit, and adjusts the axial velocity of the 
mean streamline flow at these locations to satisfy overall continuity, 

6. Subroutine SLINE checks streamline continuity for the stator 
exit and the rotor exit,and determines new streamline radii to satisfy 
Streamline continuity. 

7. Subroutine ROTOR1 converts the absolute flow properties 
ahead of the rotor to relative flow properties. 

8. Subroutine ROTOR2 determines the axial velocity ratios at 
Station 2 that satisfy the equation of motion. 

The main program computes the mass-flow-weighted value of total 
inlet to static discharge pressure ratio,and adjusts the Mach number 
of the flow ahead of the stator to obtain the pressure ratio speci- 
fied. The overall turbine performance values of ih A M,»H.P.,¥*, and 
ki are also computed in the main program. 

Computer output representing the solution for each set of condi- 
tions is printed on two pages. The flow properties at the stator 


exit plane are printed on one page. ' A sample of this output can 


114 


be seen in Table V. The second page contains the flow properties 
for the rotor as well as the overall turbine performance values. 
A sample of the second page of output is in Table VI. Most of the 
output quantities are self-explanatory. However, the following 
symbols may not be readily recognized and are defined as: 

SLINE - streamline 

ZETAP - one-half the total loss coefficient, used for calcula- 

ting Ee 
A - throat opening of the blade channel 
W-FRAC - fraction of the flowrate w between the hub and the 
indicated streamline 

BO - rotor blade inlet angle 

INCID - rotor flow incidence angle 

PSIR - relative flow velocity coefficient Le 


woth 
In addition to the main output there are additional WRITE 


Statements in the program. The additional output allows the user 

to follow more closely the intermediate steps of the solution pro- 
cess. This output is also helpful when program changes are made 
which require debugging. The extra output is not desirable for pro- 
duction runs because of the large amounts of computer running time 
and printout that result, When the additional output is wanted, l 
is read into the computer on input card number 21 for the indicator 
IND. If 0 is read into the computer for IND, only the main output 
representing the final solution will be printed, 


The statements of the program are listed on the following pages. 


LS 


CS 
om o o 


mC) mm @ pay 
Om Ome e mm 
mee “Cm eo a MmOAIQO & 
—N ~HOm 0 ON~ nem 
Cl tev entO a Netw 
Ob Of DC wedt o — ar ALE OF 
Them LOA ew em YY BO» eCiw 
mM Oca Wero I YN am ef 
Us oes EeTeONN ye nom Oo Nl 
Nem eC) © om eW —ONOd 
CO mm emf mA ON a Ap 


medtNOm OO nm Oe Aegell 
Ow OH Ora tom om ew (AJ D- (CO 
od pd come tee el cee ee [lee wOOedtatZ em 
KIO Nw etc eg KANO ole O 
COU te i OO wees te OY OO 
Saw OD 8 Mr OMeKO & Om 
<A eo eens eI JON eRe A 
Wo eS em 2c O~e WY el MON Ed 
Nm ©OemGetN™N Nem Bree ninwaw 
COm— HO D> KFOBDMNAYw~SI oe 
a—HOwdtew Oe) Be ew ODP oS OCU 
Cw atetew Cy ee om LOL 
—NeHH> OO OCNOOWS « ede 
we KAS Eee HO eS HOM ae = OO 
CY om © Oe Oe we Yow we OF OOO em 
CaN nOnrx Se CYO-ONnNO~— 
eD FO FOORFCO Ff estOoNnwwe 
COT ee St tw st EM WE NS MO 
= OO we Ore NS OO POOW KY COP WN 
TOMB HHO! eam edt eQdtes ery 
Oe KO OY me Ow em oo > Imus 
QO OWN eH IA NONMYC & eCcC dt o 
omens (Yom pete oC) eee ee] et C= en fl OO Om 
Om COW Wer ete wr lL NE eR NOOO wr HO 


a ee eee 


R1LvgALSV1, BlyR2,A7 92,82, 


Om Owe HN SAY MACY NH eS eNOS N -~_~- = 
ONS ewe ey FOCCOrF ZFCO>OF- ~ eo ee SF WANA WW WY 
(NO Cl em ete fee Yomi) eo ew & eZ lL = eee em ro eeeee e& o Co 
TC OR OKT RH OO Om we ww ww NY Se mm eA) iT] HW LOL UN aH AOHMOeO - 
OD WH Pe emt ONIN eOOOORW eo — ad il Ye | | es i I 
O> Dew OCI QO w OM SNR AIAINCA eam ~ ee keke fo ee 
mS 8 Om AA tem 8 NS OO ww ww = © ~~ mt Wim cs e@ eee e & e C 
ON <I OO w FwAie eA cn Cmnodn = won Ds tae fected Da ~ oh ee em ome = — 
AXNO es el embh Om © ew de> hie ol a ce ce Ne eel el oe ee et 
OY Kt ee = Te TS Oe eet om) O&O ow ~~ = =e ee Ye weed ee tae ee Cee tees (_) ee id 
Tee OHO 00 eX~w~w MOON eee ea e <_ Xt et et [ =) KAMAL eLL be 
RWS a> ea SRK 8 Oe HH IO HK OOM | SU ld a a © ke Alo ie de dia arom pee MSO) 
WN eter Om QO ob Tee NOONAN JE | a WONMWee dadddddod © 
NOU em Ste NI tm et LIN et Ue NI were = TL OYA CANT CeKR RR e&K Ere OKO 
Otdmtrmpwett eT MHAORrW eevee” eo Oa e Om eI NAIEINIAIIO e ; 
Ze Dee ocf ems ee Wm F&F el) eo we Km LL TO eh ee ee ee oO 


mt esta eM NRZ & CO ZNSKWROZ O NWA Het eiNO O e e m me E NIL A 
© >OFA Oo ee OCCHKNO O00 fa TN el © ewe OK ELNMMNMOMNMM | aM} o& 
OC Ost at ICO wet eA mine OM OP TT CN LNIS ALAIN ow €e eee Cin hwif\ 
DOAN stem ww Ow O me NSS HNO wm we wet DHONI OD me am ree me see mee ee em LV rm LEV UN LY EVN U0) me mee ome tm | 
HM Ze OM RK NMI ON Zw we YLZNNOOHKHO 0 oI II Om wee ef we ee _ 
SPW AKO Xe ee I JUS NE ENNO ANGE QOOSCOFOFOOOOO00RSH! 
Te Ee DHOOM NEAR TOU WUE NW Taadeaddddiddegdyaqadqadaqdcdqgd 
UL Li el Oe et eNO IO IF MO OK RWOUW WOW OW a 
SNDOFODO~CA OWKNOFO> MOQOOOWW ELEY KE XU KLM Yee eee 
ABAORANMALNOFDORNM BNM SIA . ~ 
— mM MK 


12 @) 
~ 


t AGN@TO-ND MOR OTOMNN THOM EOOM 
oO CQ OCOCO Sel Had at IRINA NANNING 
eS DACCODOCOCOADCOCQODOCCOOCOGCOCOCOCOCOCaO 
.O DCAQSOEOOGOAGOGSOHSAIOOEOOHoag 


; 
' 


0001 
0002 
0003 


en er ee te a ci, ee tn 














116 


e - mm 
— — — wee 
me & Mae 
YA YN ae 
@ ~ + 
oO © —_—om E 
~~ aw WU) 
a coal aaetiaaeed 
e = 
a om mm 
A NM +4 
e e Mmm 
4a CL eH 
~-Ch e aA CS me 
© @ oom >» > ww 
OreG o aii VOY 
oe eM UW Ww we 
m~coch e QA Aa ## 
ee ef) eo _—-_ 
COO eo a NN tt 
ee en > > — 
LAIN UR - J 
ef ee be ok =—lr 
CSC eOO Oo O +4 
eC) ee ss = WN 
— ein - ~ t+ 
aAncc asa WN 
OatOwo <t {tf 
*L) ee o e foe n aS) 
=“ eA) a A ww 
<N 1M Cw wo Ne 
— caecn e ~— we 
a em ee << <t HH 
eC a est =: FF mem 
UN ot aIvucg Cl Om 
~ -_ ~o qacd<c on ME ce dined 
C: a) — ect eo © A WM wet) 
=a ee mt mM eon o e OM 
am CO ee MO ete em ee er © + + 
med tN et ee GQ) eAlst e  -—— 
Veet, an mem ef ee aomm 
mt Oe DO mie et ecm & _ —— 
ett wl ew th wee em ee Sto “Max 
meme Clee © FS oP —— * ~- wY 
—eemee PLM ASA MM + sts KH 
ae EN pa Si? SO ON ke <T A UL <I p= b= on eee a  _—m 
See Ae COU wu Wo ao NN 
awe XS NONC ZIM ad H a 
wer et ed wee AJOL eee wee Cm ee ee Ge ae card Om ee ot CN! 
ead mee we MIN me eer AJOIOLM 
STP meee ree fee wee FT OOM CO) wt tt ee ee ee OHO + + B 
Cm mOMRMON DRY IS Ise ye KC KNNN ean mnm ne Qi qadqa-=— anna 
=A MOMs ON CYAN NALY LAM LAMNHAIMNNG «GG SILL LL OLLI TA 
eee ee ome © ONION eee ee YY Hee NN 1 www we 
LBL UN UD LT be ee LAU ee O00. WW NOFNOPNMAN AY Yeh SR BINNS HNANC 
——Seee TS CmTOCOONODOCALLYAE OHO COMI We I «KK «KDE dK 
Pil Neale (el ay(a ayaa Ce ee et] i (Yf] a vue nope en uu Sei oeuvwwu ud 
teed ddd i we Se CO COO RMN AO QOVOWLUE LLL Le miner 


MW WWW WOW ee Ke KT LOH YLAXHNYNNAN AES Ea IHN YNNNYCAaAMNNMNaxvaia 
LOLOL CK KOWCURK KAN NENA CVC KCYQOUUTAAY AAA oY OK mK OK OK OK OK 


—_ ex 
st ies) 
fe © ~~ 
~~ 


TOO ATORNM TNO 
ODODODORRRRRARE 
COOCADCOOCCOCOD00 
COOCACCOCCOCCGO: 


RM FO OMDROGAIMEFN OY OF CHANMSPINOMH ODOR 
MA OOMMOMM Ss TESTS TT UMMNMNWIMNWNDWN OO GO 
CDDOOVVOVDOOCOCDOCOCOCCOVCOOCADGDOOOCOC0OGd 
ADODOOCOCOCODDCCOOOCTCOOOCOAGDCOOCCOCORG 


0063 


117 


2 ne ete SS eee aaye ee - 


WGT pV 39DELY39R3—9SHR_, TIN, 
TAR 39 WGT pY4eDELY4,B4,SB,T10,% 


ry 

> os 

Cc - 

wt =) 

we o 

- © 

oO -_ 

YW - 

Ld oe 

ww @ 

oc ww 

to e 

mw oO 

> a 

al] eo 

w wWO 

Oo > 

- oJ 
-~ wm iW | 
Ws > A 
NM eo e 
~< tr Co 
° Oo > 
oa) x o 
Ye - _ 
x et Ww 5 
eB = 
—a<_ <q o 
Sh oe ac 
<u Wl Ld a. 
eA La La] <_ o | 
Oe Cs - -_ = 
SOU 1) A a | 
Sats Nae ea . [o¢ 
ee ee | o C! NY A) fan] 
Na wo x ow a it le 4 # 
we - o e = + + * * 
car at aq ¢ << Y (A 2 ac 
ai f= =F FF o =. | = = i = 
ce a) Ge) Ss < <q < <_ | 
OY tet ed be 5g * + * j 
Or HMA WA WN ee oat AL x las) 
lee - o e ~ <q oO << @ (oa) 
ie; (2) j) (& © UO Looe W WW O \ 
~o - o o = O+ N +O? + { 
ae O.5O “© a. MINY | * SSe fend | 
- e Cs rs @ CN es oO ‘xX eZ = | 
+O Cc Oo © Ue) Nod = aod <q | 
oe e - = oT & - Nm + OS ~ +t ~ | 
OsvyeaeTetfeyr MW Ned Q eo en Net OT WN an) | 
TI eOC 0O eO & fH fF & rod Ow aOm — @ — : 
rFOROAOHeO Ww WUes fa Taavea) oxy mora z=. @w WY 
WAdatdedew YW ORT - + Na + + & + fed 
‘Nh ow ow ow ZF Ze eo oe wh News NNO TINGS °o 
wt St Jt et IO e © CWI w ema | ANO W DOOO = ey WwW 
COrROF- OKO ern Heli N~a J o<f (i O tJ) Smee Wom Mm 
AACAOAVUAACTYW IO —~et Ww mn nh eth ouncnoool~ |\eOoO uy 
eT CO +O 6 ow Oe ow CO-0090 | =O-0-07024N4NNH=MNOO ot ee OOM, 
IYONnNONONAw W == | MOMOMNMOM | MAMAMEY YY AULULCO I amet ec ec 
md hd cd at ee be CO ee bm +N we Bw ww we KH NONNM © omit eon tt 
awd eye) el eOTTOPOOZ 2a Ike Cee eddie Ie IM SMe HN KIM 
IIs dt It IOC <I il tl Hwee ih th | Ue ee we Ln et ed wee re ey ee ee 

OILWOWOWHMNNVOHNNNHARCNNNMN 


OUARAD AO ORI EEL soa eee ae 


OQUNONQNOOCLOLUZZZmm {OO gq@ Om Toda DQndownmaxrxexOxaCuu>>>>>> 








om oy = = 
+ nw oO fole) o C00 000 0 OO ro) 
Ww am ££ AE DOK © BO ran 
~ ~ | NOW “a uw“ 
, | 
Ot O BN MENOKFOROANMAGCIOR © Onna ENON OTOMUMLAORDAO~! 
H~ DH D D DOWMOMDODHRAKRARKAGAAABROOCSOOCO0COOO Oe at at at ttt HO 
CO GF BD QO DOQQQQ00QQ0OQO OOO ea era kl hd eh at at ald at at peed a ed oan oad ld pad ed wd add ot md 
ono oS 


. 


DOODODQOOOCQCOCOQ0C0C000 a et ae. CO00C00 


118 





o e 

td ~ aw fa 

o uw ~ 2 my 

LL - NG oe > ae Eee 

WY We a an Na ue 

~~ 7) e x Y) Ne aw 

—) le <a UL uw <<fe <I 

NS e Y mt Or xX ~ 

o<_ re < ae NAO} uw 

WLLL —Y rn = Se) >a (Va) 

> (ae om eLL ON a —~- 

e< e < S16 wa Ww foal - 

Ne Fe) an x Toa) Dee ew 

mA Nien WY ow >WV ae ao nw 

“NJ om e) —(CO oa J LiL o 

t Hn uc o i o MN eW wm Wt 
asaz WW! = aN) Zwo Yam 
mA) a <a (Ta) =~ oh © e aaj 
NI o<y - am OO. UL uN <a 
rm < Oe D4 ~ oe —™ CY i Mu 
— Oe hk Ee NH xq SU 4 a em 
<W — . ui < 3 op © be << 
i ><N alee as re 8) oN] ae <I ao 
c= odo LU b= QD UW) ae <a eof — We 
ay = ams ee = He >xKo ta. 
u > AO a) N A wi oy) & coe cell oon) 
‘ay e<cl oa ad == a Gi e&N za bol ead 
< AON Oe ns LO xa mm Ce 
= Q ec bo ond wa) tk «XY [(Saiae Cra CY 
a eae ac BC bt AL & Wee a aw 
< a<J lL! ~ o oS a QAIN o ~ oe 
oe haOAy us Cc Li, em “cr Ws ibn 
ci ene rm or w ~ Tw. ah me eo 
2 Awth A: = LL ~ van =§ MN ee <a Q o& 
e =< oo -4 aA -_ a oh Ons Pd Ss 
> ante ~_— & = _ Ion “Na co — ao 
x ee eo bow bij o — bmw rm ONO S-: mt 
ts Cum za “ 4 Wot aee Yr B= 
4 be UW) Cs ae) (oe be OW be > poet med Ww ~e 
Acme Cok Wu ~~ emu) Aam ake AW 
ie eO & x oO b= =e thi «CO Wi «OQ 
ia’ C es eae FF SC fad xX & Is « eQ co a) 
= HKnn MAO _ a je] = Ow Owens OFZ WO 
Ic rena Qu:co -«m @ ” AWM He~O ale Ae 
hood eee aae © e + eA aE(nd Lis wal Mey <n 
a aN rszN~ FE O ro eO oe IM e eo ow 
e >< mt QE UW eG OO CC Oo a eM ame ' Orn Wu 
Te: Ho NO om (4) ° aver NEI NRO SIO 
= pe Y) egqt e Hee Qa >vw qa Mey ez 
< I «Oo —a>—-4 «moO re eee hm LL Ike Ne 
~ UL at h— eG ~ es net et ot CY et WSs OmO el a pr OY 
ic ah be azo Neo Nao Dawn ONY Naw ae 
j- ve Aeq ade OW am amHY 372rAae Tr er sn) “a> 
j—- oN QTe ~iuj @ Nw me eee elle QA e& 
| mae era ee ~()O Bore CO Qn WO CNS ONS OWN wZ 
| sO evHOOe EN LON NO Seete aANM AeKYOramm aad 


OL Ih hm) fl beh Ss ew ] 


ZzU;eCocoeco *# OOLWOR>~—~ Cee 


HOeey! enozr! 


<S MMID HAF OM We ole ee eee HOU Nee MRM eNO~OOW 
oat e>y IN CDUNON ARR ae ALO eweDOAMe OCOCOMO~s UZ 


[aan 


IO OO I EO eee eh ee I PU 


i ee 


l 

‘] 

| ro) 
oad 

| cc 

| 

{ 

{ 


L) 

Are 
TIN TAY 
el | ol 
O000 


a eN eS NO LON mt mem mem sti KY mmm tiem UP HO Wore OO 


Od Wt & 


at od 


AN 
=~ 


wor aon 
NANCY 
ot md eed pod 


Od00 


a oe UC ae er a ee | > A 
<“<"iCKINWCRrRAYULOMTNLS YMA YCOUUCAYNACOZNOAMMOOdAeNaUNeCOodteg 
CZCOUDOOAkKONSBOULE > >>> > eFrUL Ie eK O IMO YDOOCU>F YS EUOOD! 


Conn 
Men en 
mde 


ol @ Tas) 


QOeN 
(a) ey 
© ad 
cc 


MHTFOOM DHOANGATNOM ODO 
MOAOOAAA OST PP rrr sun 
ed dt de ted dd eed td ed ed ed ed 


OOOO VOO COCO COOO0O0 000 





wo. 


119 


ee] 


© 
ia) 


OOo 


= AN eCOMKKLLOZNNCK INN eOrvae I it mNZee eHOekeZYSZS &Ke Ue 
! 


AION ST teeth INI COWU NN IN & Beeb Se “| 


! 
I 
' 
4 


ANN od 

Oo Oo 

NO 

fan} (aa) 

WO ho & 

UD LAU uA a 

md cand mt tall od 
'OOoooo0O 


! > - 
| o ~~ Cw [a 
eS uu N ex La 
UL o —Y oe = - @ o 5 
WY u ec et A Cate & LL © 
Cred va Ca KW Ne S a ° 
—- @ 1g LL ou MI o <I lo @} 
As & Ta) —a & ory o ~ 
od ve <q Ze NAW u.. ~ 
a WALL WY me I7e) > ee WY eS So ~ 
_ > [ag = ou aN fod = Va) ' 
4 oat - <> me NAW. e@ « ~ | 
Oo Ne —Y OQ zea eo —~— © fon) | 
” aN wa | «|O ow Sue ON w 
| c diem eH OF mae oe ° ran) < 
rs en Ue <f o j= ©& GN ew Le (a ae 
= a ee uw = on zwWo Yam + a 
WY mA, & b= @e a WY Cr o-— & <a Y + ‘ 
i ~ Nl x oct o am Acai un <I ~ ' 
’ «x em <f Oe < ee “ Ory Mu x< -~ ' 
& ot ob Mm WO Q New Su. oO Cc t 
uu Iw Ka LL <r oe & be <f Uw ad 
ag >In ory ae Om Se 0 ee oe om 
eS o<t > UL! a) Us om <q ed Wo Cc - , 
- — oo ke o Dm >axew aoe e < t 
lan >t el Ol CJ Aw oY) & eee fan] (od 
= ade -ox > =z N eX oh = a 
o AWN Cie o LLC oe «<< OD Oe a seo H 
0) QA oe<r b— om UW! hm of adac amae4 uJ a ji 
ae Opel ral, {a ba may & Nee Aw pa on) j 
ond TLL ee o = — Qasas ee co Some | 
oO hon Wis en) LL ote aa WIN ce oa) ' 
ta eee — co we =o acl ke CO fe @ fee | 
uu Mm fo J) Us oY) <I x eo ae vy) < 
(eo, =<an ~z al, “ob Coton = a (a d ‘ 
= ° aaev se =z qa New a <0 = rote H 
o aN] eee mL Ca =< & ANWO mt ond ad +> 
oN ~ Own zQ ml WoO eee z= N -~ 
(e7 -~ ee ef) -C a Oa” ke DC pnt emt eo ‘od aN) 1 
Liem os Ada aa) uJ eOW Oc ca Nw Wem wed ' 
SC AN = reO¢ <e a OD «OD WwW eo xO an oO Oo 
on zero ~w~ C ee eas EF DO od wp TY eo eC X%aot m= AN : 
aN] eeN HAWN acc ee MTOM OC ott om fo i he) oe oO aq o 
ee KN | & | eh ole 4 Awc x AWM ea co a eA& ACO ~ hal 
C OeOM + eee Mae OF wmaere Weq WIND Occo aA © 
Tc i SwN = any ee FF Oegere. THe Cray = FW ~ SN 
aA Accn — >< oes CO UW eMN CE Des oM ow ST Ol) —CO ON % [o \ 
o OO - w= om & AJ et CO em ecf Naw ena oO ON = em 
QO Ww Cort Om mp) eqde ewnOrordt Gee OeBD WO =ANO | 
C K +O Om I eA s—>— OO eo Km NA ee dx er wwe eH ' 
CJ CAIN ww tL <r ef & BO SOY US ms OD mm 0) p= A OD PAI We | 
ww | al AZD New DAM Mer eM eae SCOD Ax 
~ Yur —-a Madu Xm FT eewedadt Ow PK SWC oW Cua = weal 
{| eo wr & OLf\ he | wet eer CY LOY axe ~— UU Ome & wl edt) efQ em ww © INO ) Be 
Ics HAO ewe eo eo ~OO QOOCOmWwO eA OMNO~Z oe HAO ot Oe 
NS WOT tO ere OCOr EN WOTAND Hertovrnweato wWworraqdano 
1 au} Wow Nee = ej a O3«xe [Or eg~eZzIN Susi ta~w © 
Wile elles HHO Sef D oO Oe OWN YR WN ee | OU} I Wem elem Qnm 
iCDdAWLOo =—a@ eo lo CODINOOMM ¢ OFCOMnNYIswDYILO &~DO - 


NZNO m= | HN MOOR KLLOZNMOCKE PHOekKE ZOU eOZHVORK~WI CO 

Om COM OC Belem NOM BO WHO Wee Z WO LONMWN I OOO! 
SF SON KO pe IID IN) © om SO ad Se uo 
1) ze HUE wail wdOwd Wu e Stee =e QO od ON Swe ae SZ I< weet ot! 
luGazaucuanoanlcoHr<cuoa <x OCOdaNCUdxeuUOQOdwuUOotreusOes 
(SOG E=BOGZOUGRSaAr ONS OOU Oe BOOUSh emOx OMe MOR AON ZOZ 


a mA wal aN tA = i= 
{ (Sa) oOo = © oO tN Comal Ve) O Ph SO fr Oo © 
an N WN @ We) CON! Ne) a fon Oo AC a A 
ot SN WN co COCO Cc foe) 0 0) te off 0) NN oO Oo 
e 
{ 
‘Omn OSTINOM OD OCOFNM TNO FD Ox NA FTNCG POrtonnanswn 
iO0OD0O0 OSCWODOOO acl eed ined ee Ded Sel Doe Doel Ded PMO Cm Dowd iad adiped ol 
pod pad gamed aed el wel od ead weed oneal pam} ol coed mel ld el pl ool el owe) oe! = ew mel mel meet =) A ed =A ed mel 
‘Ooo O2O0000°0 oogo coo |co oOce | a ibis is 
| 
H | 


120 


{ 
' 
aad { 
— WO _ a | 
(as <I = = a i 
a a <I = <I 5 
ve a = — | 
G = o La | -~| 
¢ =~ ome =| 
— 4 — 
= booed — (ad 
-~ fal, — = = -_ @w a 
re (ee <I = w 4 mt < <q <q 
© > WwW =< tw z= ~ Oo Land = 
to fed a _ et lw Ly 
lead <q = VY b= La} La | 
=) e o =z ae | 

al - cot <I a=! — 
ewe Ym — i e (on) = 
OM k& qa ~ <I —_ -~ ” WwW — 
fom aq "N = es) —_ a - or 
oc - > lw ° =a = (ed a <i 
c: Ww) o (oa ~ = = OQ bk 
<a i — roa oO Ls 
' cone — = -~ w | 
ke JT a <q <I om = 
Ww I ui = met -_ mt 
NS UJ <a o ty > & (= oh 
= Ls) e CC aal 
ern o -~ = Le | 
pie e ae ~-~ — Ww = —~_ om 
-- UL — — — -~ = « m= CU ~~ ow 
Lu - «~ = a = “a NS > ~< we OO LL 
1) Ge < = <_< ee) o —Uu ke} 
wr e - - ueop7z> = fad N ek! 
(ol (oe LL - jj} Wwe — a hi 
ou Not = < ~ — C — om 
} >} ~e == e GS a > b= 6 ee 
at eo oC) ct 4 oO (Oe =O > cen QINN 
_ nD a ~~ oe YW? ( ¢ em iW = e - ¢ WY ox< 
- us Com” + errs NN SG mm oO HN eo 
~~ oC ~— * Me tan LAU OL > et = — om 
Oo N |W mN Cale @ Cm UNE Hie Nw =e 6 p= COO 
at OM A ~e W =u UW ow 3 om m e = e* WwW om 
~ « CC Nes Crete mittee MO ss ST [Oy ms Aa 
oe ot ef eS a Sheer UL DZ >See = | = tft ea 
NONON Qe ww WROD oJ wee CI ge mad od od om ox Q ) J wwe oe 
at ON ON on YW we WN NOewelr ee NN J me eo e WY ee ee 

i f-™~ en om —) b= *@ —OOUr & wn oO-— ond NN Oo 
_ -LOoLoo— em Tit etbe DT Wen ZW Om NOLO mL ml 
§  DULUINL A Detem Tm Te Kt Oe NMO ow Creme Mm LWINT ONO Te-oc> ee 
ee OT OMOE OCOD MT AHN CFG ett ii i Cat etetedt eam OMO tN ee 


(OT OmO™ om oT OW om ON HT ON HD poms oT MINH oT mw oT eee IT SS 
NXE | MHF Om Om ON Gm Om ee Om HP mmr ON ROMO HK ODDS SI HHO 
ISFOZTO RR OO Oe OO eerie ok OO = — Onenp- 
CaeadzaZwadwatwdwdaowdWwdgdteenase dwaiwnwawwdiwwwi|g ewaw 
HON be Te TR Te Pe TNE Te DNNNode De BMH DSR Re PR eHR—RTE HANTS eH 
ad Oh 1 lle dole Male Saale M10 a eae Ae 2 - Oe ee le A a Alo Ml elle A ee? is A 
|jQu Fe Sha CxOrOrC x 0OeOxOrOlamswae-Ga0oa0anrwOreeO Oxr0u 
Omelem qe SU EI EFU BU SULT EU SULAO MACAO Feel BPUOZU FPEUTEFU OUST 


: o- om O- 

1 (QUN) sp GaiePp fA Ip ee Co ie] ia oer NN OO ~ [oe] ano 

| asa at AO “a On ~~ A as Ne = ro = 
~~ ~~ NO Vea a) uw oe 


oe 


SHOR GCANMPNORCTOANMTNO MOROFBNM SINOKCHOAINMGATNHO KOKO 
COP TDCOCOCCOOC OO CHa ett SRRINNINN NNNNINNIMMGMMMMM Mmm 
AAS ANNNNNANNINANNNNINNIN NANI | INNING NINN 
DOOVOCOOCOOCCCOOD00C0O000 BO00000 | COOOGOOdO0000 Oo000 





121 


o 
tad ae 
© = +4 
N ro] ~ N 
! = = - 
oo — | 
aa = | 
—— <q me 
i — po Cael 
‘ ~N N Vp) _ 
H < =s) ~ — Y) 
: b= - N = x- — 
: foal — a ri LW om | 
' ~ (ay OQ 
i ™ - + | 
t N Bs a 
<I - N ~ ~~ N 
LL Ww => -_ T; kX i 
—_ (ec as 4 ome = a 
<q ia Sed UJ 
~— CJ (a) i 
— mo oo = % 
N aa ~ oP, aN — 
a. et ed i) om = feu 
~N Oca fo} oe % om et HE — be 
eS tou. > law © = ee om — land 
- ~ J = ~ pam) te N 
oe o~ —~ Qe Na Ww 
lies ee ee - oe = 
~ A — — — Wem ton 1 - 
a. am > NON We #N~Rr ~ - | 
e- > > D> FANS RNR ~— a [ok } 
-_ Oct o > + «Ome —_ 
- wom ~~ % om LL Se cee fe ceed LL SA tee “AJ 
— ome Oo. Cl am | an MH Let | hh ee fad 
=... mo ou mee HH Of ee uw 
Np oe cn — an te LW ee Qe ae 
fa) Tw ANC fa ue =) N Nee e ere NSO zoQ — 
— > ce elL— = o> D> wm OAK Hr =i at ow 
NIE ME ON ~ { Q. New k ON wee em et be ee HIN = a) 
| mt SH HN © on Poet Sa at fm SA pele me come ee hm Eee > = ° <« Wi ; 
Ww wef ma Oem Li CO me hee mew ome A) | O. — +m © Oo w i 
oe 7 eT ae ce De RR HANS | me COBuUn wt tm | 
— SaONDAI O&Nw Um ~z > TDR eR New ee eee OE a—een OF- Fe Ww NN! 
Va) ~OUMANEA YO wa HH Jee OH NL ~—~7xClOre#k wd ° 
m—_ a= © ect — eo mn LL He ee mm ff ere i we LU™ fe ¢d tre 7 ioe 


mm TN Tm Le Oe mem HM OOW HO FS CH Om TR OWE, KR IL: 
wt OO Oh mer TOL mUIa SS oll owe | me ebb bb HO 8 em Lm Oe tN Tt et oo 0) 


~—NO sete ft ered NOe ate Qh ADOw wR WR RRS (OOF OOWMFN TIBETAN MHed eo! 
Neen Wy NernmtN ew eh FEN Ihe RNR Rew NY It HY WeAMAcwOr~rOndc enti: 
Ot te Det eee HE AO HE HW OH Ie ee me OR EEO me 
> Pot ee ee ee Ot —O— =| Pera mem KO oranandax~ ~~! 
e RPUO INN Rew Fe OF we Omen ee mete Seed O OST HN QH ae FOF: 


TN IN I LN LE i LL Oe et et et nee tt meee te ee ee ee Fe LL HY I Tmt ] SUTIN <T 
mt be SF OY et Lee el FR YE SIN SRE ON Dew Nem YLINYVIN $+ LOR RAS eE SVS 
Oe eK CIN NN ce oe OW IO ETORMaA 
YOO | OCINnwIWOWaA COrFOCA=OWR RE ANEEEOY~NWO WWatrwotaucwOoo 
>FUBCA-aAOvEeFFSUNOU FWCOOrKrAAkKeWWWEYdoOOUOLId-OlLde FRUOW 
tae) = eat gon] =! ' 


fon) i uw oO (o>) © ~~ © 
ia NO NI N= oN, wf ' NON 
-_ N NI fa] ; 


SN ASINOE OD Oe NMOS NOM DROHNMSNHOPDHOANMIIN OM OROKAM 
tr yer ye NW NOW NNMNHNHOQOODDVOOOORPEERPPEE EEL ODDO 
NN ANNNNANN ONS ONAN AQINININNNNAINNNNN NANNIAINNAIN AAAS 
O20 DQOCCOO0O O00 OO BGOQWVVCOQCCDOCCOCOCOG0O0O COOG009 Coo 


| 





122 


U. ac 
Li 
7a} od jag 
—_ — 
w = 
<x = a e 
” z Zc wn 
mes Ww ~~ 
~ x = es me 
i << oa) eet =~ <I 
tor Oe ia = ! Cam 
1 q = — ef 
IN & tomes U u. Om 
ie WY) = Wu oN me 
je «2 Yr « Naat oO 
|= <I Us aoe fe - 
' -_ mm 
a Ta) Ww m~FPre © 
YM & in a YN r—_— + 
— tm & a x ae Nm fe 
nN MN -~ ) Cot et 1 hot mer ~ 0 
- & = ae x= 2 ectyw - 
e ~— CG etordyry on 
-~ fad ve QQ al + ] 
— — Wi Fu tae FT 
— ynWY a ee em (NY) 
(m  (eF ef ood ond (1) wer CY ~ 
ee << ce =z om NMe NV Te) 
QO bf U.- Ww CG Amare OC Ur e 
~ tt ems fe ub Yea«yn nN wo 
-~ CO < sr Wee te C te 
ont ew a “~N es elLey -& iL 
Sad cat (eS = ZaoO+—aI A We 
~™ LL. a c Combe fo an 
—- WN eli ve aH FSainwW ee 
—- << ue hvala a WoO 
~ - - LL a ee ee VY) = 
-~ Lid fo OF) — eI y ome LL 
—_ LL = - O-ibk-<{ e we 
=e oY) =z Hm FT etl~mt ff <tin 
Laamel td ew) ood eOwt + oe 
— tot o-— =: t zat Om WwW Cas 
a << COW S me Ou eNN =~ 
- - we ee Wo ad ON + ww Ii 
com Ne OFX oe 
hs ora e emingd UWS 
_— ta 2 ON me eOO 
a om WJ eC OCOurnnn be Taal eo] 
NY ~ << I= lS KAI 
w 2 10 Oe WW £ ANAY + Au. e @ 
QO = elu! Ce Cm & om fm e ©«OO 
~~ x re ret AWOQAOW $+ = wrong 
— =~ em ~w Wo DLL oe I ee 
oO © WwW) Oo HHO OO ee 


—~——T INL = at GE oes i Wu Oo 
[o@l Onl Ca) eo NM TOS MO oe NENT mm ec ON 
SION APO OME OM OTDM Ee SN (UO 


ee oe™ i ewe Uj eld Yam af Oem 
COC eee OC Ome OSE Te ON Otte 


; a) ee e eK HD ee eros 

WW ODW WUE We UTI NeW Zee 

Re EUINER XRT RST TRE SONI TYE EH Ze 

vee OL ee ee Oe Te be et ee ZO 

XXOOCKC CWO CWOCKOCrYNCOrFCOULOZ 

SEL Ou®BOFuU FU Fuwdty Fun Few ow 
ped 


od oo ~~ Cre 
ie OS mO N tT ing) F's) 0 OO 
wm a mM O Mm (59 aa) oO an MO 
= > nt 


TNO FORT OH NM FMWOrORO aNuMTNO 
COD DOD COP AHR FORAASAHAS OOOOCCO 
NAN NAN NN NN NNNNNNE MM MMGG MS 
Coo C00 O00 S00O OCO00090O CcOO00o0 


ae _ ws eae te 


viru + Ae Sao, VAC 1, VAR. VAP3 I 
2 
A 


i) >t 
~ 


ARAB 


Ow CHEN 
> HKHEMAL> 
Ww OANA oe 
222eACeXIw Ht 
SOc gqoen 
p= LPo>>te- 
DW Hn HA 
CZONNNY~aA 
SMWwTANNSIS> 

Cot EXAXa> il 
DeOddd lice 

NOUS >2>06 
Oe 


=NOTIOAOM & 
SGOGOGO000 
DOCCSOOO 
COOSOOCO 


mx 

CL a 

<a> 

>t 

tN 

fee 
[ik3 (a eieuira 
a to 4 
=> Catfeakts) =) 
rucoorc 
mutt tl wl 
(any an X= fe al Gb Yo a UE) 


DOINGS 
ot ct ot td od od 
Ooco0ogoo0o 
Ooo0go0o°o 


124 


ecoa09 GOB OCOCO0CIG9090900000000900006G0000000090000909 

SNNNNY NNNWNNNWWNNMNNNMNNNHNNNKLNC“UVNWNNNNNNNNNNMNNNMNMN 

a a PS nS Oe SS SS PSG OS DS I SO De DW Ow eS ew ew Dee ew ee oe ee ee Pe ees 
o 


(2) 

e a 

—_- A £ 

(Va) > 4 om 

os oe om ee 

—_ =F Cw 

-* A OW 

Oo Kx wetY 

"nN A ~—~e 

~ x Om 

fee) ~ coped 

- & = ww 

~ «K OW 

Hi fe Oe 

tL! x< = om 

Cs ~~ wer ed 

~ QA Lea 

~- «Kx ae 

ee & O gan 

ie ~~ me 

~*~ kL Cw 

N «K OO 

uw a om oe 

to o~ ee 

< (og LL om 

~ Oo 

a“ DN —™ oe 

Py >< ol 1 

[Ges O eed Se oe) j 

— AUOw 

YW 2>nT H 

- Cis eM ' 

Cc LO eo ! 

pall) eo a= 4 

ee CE awe 

x< oc — e 

Co Zan> oo 

bot coq e eae ond ' 

Co o> = ~ i 

=F A omni N 

— Ope ee ee — 

-~ UWO Nb NS <q 

2 eel o@ * — + ' 

yxy IA om % nt ies | 

e *0A =O oO ~_——™ — Gs ! 

= Ox an © Some bea = ~ ] 

— Oo — & © "ac! ee a be ! 

— FPaAxKO va > i § t—— + ; 

Cn. ee en cn aa AR peed Gee ome OL | 

QA Tuiew (o>) + % OL we oe oN O 

oO x<od ~ ~— Tua Aor - | 

Vi YOM ea ° Cc —~— = tet tee Mt + Men Wie 
SN Net aHO os wt er Le | mt mt qqyzyeo; 

a wt we Owe O (BOT oO we NA+ OOO wow SF #HYW &O 

Zz F2EZ ew a =O Beth Sexe « Zzeaan QA —— Ot 

— WO aman OWwW CHHOO 60 Catewt dade | eH HUE MONS | 

— a ome hwy oD i | Dl Om ENN OMNMST 1 OCKUA KUT Hew |] =O} 

D OW Yt SeaeMe7Z2FZ OVP H-HOO tt ~“OMov¢s LL Km 8 oe KOKA NE +h ee NS 

Det I Z~O ) poms. pao ePHIN tS eNKOTiK—udwdHecomnaacnwen#eu wo y w 


MISO eSeS oF COW & OY Omen eMR OWL me OM tH Om Pema ane nm NK Ore 
meta tvs ome A tO ee ch rt fl mete Ip et ee ee eT tN Nt La OK md I met A et ee tt 
DOW JWR wOUGDKO00 3D ~~ AODOULNWVN ONES ee SOO Qwe KK Ow KM Kw Rw Ewe ww OU LY 
NN ADO QOVOOOOVCO4IUM A IU KOM FBTOBUUWAU KMReEAAAMAAArAMYOMOm= 


— —_ axl 


rer DD (=) Qh Sc 2) Oo Oo 
= On omd N Oo tt 
os) Oo 
= oe 


TMOOH LTC ANM TNOE DF OHMAMTINOM OPFOCHAUMNSINOEDOROHAMSTIN 
DOO DOG wi rt rd et et st tt et IOI OI OUI IUCN NI OED) COMO MEN OS SE 
DOODDOOCCOOCQIVDOOVDO COQDODOO DODO COOOODOD0ONOCOO0OOSC 
DTDODOODTCOVD DODO OO QDOOO DO CODOOODOOOQOOOCQOCCO0O0OO 


| 


| s et gee 


ocr: 
0002 
0003 


125 





MANNMNNNNNNNANLNAMNNANNNNNAMNNNNMNNNANNNNANNNNNWHNNNNANNM 
eel tN cd ced wd ced ed ced cd wed cred a md ct ce) dd ced) td od) cd etd ed od eed edd) od) oI I ed) I) ed I od ed I I I I) I 


‘ 
$ 


—_— 
eo 
a 
! 
wz 
— 
<f 
5 
<I 
om J _ 
—_ WwW 
rou 8 
N od I 
am (<< = 
a Mw bond <a 
— a o~ — 
= <—ot “4 ae oes { 
a | ' = x | 
: % Ivy 4 % | 
} <I bh om — ~ = ‘ 
| an LU em <I CO — 
LL! ‘oan | * S-9 bh 
oc iw <q * N 
I _—— ac — © O ~ 
—- — xc Qa = —~ — oJ — = 
-~ Qa o | + no ~~ Uu ea) pa 
| ~ I= InN © | x Yer el ~ Y 
ie =e, OS Gus ~~ — < a0 ~ot ON F& 
iw =—t 06. wdHtm A Ie ad ! ++ met Nee fem mm =D) aw 
} aaax Iqa~ a- ett © — = meu qo 7 9 oe) a] 
“(O, ex xs [=O x« a} uy - =x xis OW 1 a) N 
Oo ~a¢ean e+ I~ —_-— ~ we BLO N ~~ ~, 
— »5»< AN aaa aA Ow ~- ste dee OO LU ed >m <_ 
” w+ OT %#+a ¢6§N§N JIIOO ~rYrt FPWerOnHe &— CWOrOANEy 
fae O 2aax@éan OF w9UeK UAaATH oN ow on ee > HONS eNDeaw « 0 o 
™- © eO e# KAQRKO © &— m= KOe MO I me NT RHO HNOOT OU DNOIO Ae 
—_ 0D emi AK MA COHR~ AAA AAK HANS NNN CONOR me KIN HY I 
l= —OUCNWMM UH UNKCH He KOA PNA WAS | ROA me NO MO NDC IH me OU I HD et me 
J NST UM HK HK TAN HAH | NAM NTAayvY ! pOVOn ron HNNt Cot nares 


‘TOQAASRCKAATAUAOCO KLARA Hee OY YN HORM > ONMNY Or PUT ROMO 
CL CD KE KEL HK DK KE OP OK OK ee Oe Ot wel I EO ee 
‘OF AA~OAGEAUWWADOK~OANA DHE Re ww Ow Om Owe WU tO De OO 
pL < OOx< «<< TAR OOAANAAAKAA Tete YORMODNNOONMWMMHRKULOLNNNOMNO 
\ 


oO wl © OANA OS oOo O0O0OW WO oOo tN eo 
oO so — ODDO Ht NINO DO ~-oO WW UW oO 
jo —— ada eiadt = md oOo OO ff 

st ot old od 


— ere ete nae ee oe 


ee er a = 


OE DORON MINOR KE DOHAIM ST NOR DRORAMTNORDOROHNMTNOR DORON 
PTS TIMI ALAUWIAIAUNIUN OO ODO OD QD DORE EEE EEE EP DDODDDDODODOAAMA 
CODD ODOD ODOC OOO DDDO OOOO ODS OOSCOSOO DOOGVOCGDOOOCOCOGOCCSC9OOO 
OCOD ODODOVQDOOOSOD TCOODSCOOOOSOANGOOMGSOSOOCOCSODOOMOOO@GQ009g00 


| | , | 


126 


CAGOCCCOCUOOC OC OCCeCCODBCBCOEBCOCGND00000C00000000009000 
NUN NAY MNUNNANNYNY YN AG NNNNMRHnNMNMANANNMANMM AMMNNMNNNNNNNMNMNY 
a ee eS a a SS i) JV gS SiS i Oi 


~ 

3 

+ (@) 

~ oo OAR © 
= (Ep QANODOG coonad 
” CG Sc SeoQaQgvortroace 
x COCA e eOCO eCO eQ 
rd Cc) Cc COO CLM ee eMOCwW OO 
Y OCD CQMDOMO eONINDDOMAOWM e 
+ () (as) COCCOCLE eN eLOeAnin «VV. WO 
~ GTS Oo QO 7CODO COC HHODNONNHAONONNWN 
= N Oa OC CO eO eO eS e eNXRNH eCOCOnnanncoron 
em WL oO e OO COWVM DO eCHODDO BOHN DOQOOHOOCQ™ 
AML © BG Ge DW ett OQet QW aeH~TOrnenccececrno occ 
ee pet be @] ~ ee AMON ON AMON N OCOD AM 6 COOO @tNnNO OD 
Iw o- ~ CC DN OCONCTOONOCOOOCL eDOnMO o pQGTS «OC 
me CSUN = oy 


Q SO ACQOVOQVGOOCADONDS *ONTDOUGUMOM « 
+ NQWOOOGODC OOTCAC COOLGOORGS *OOG °GO 0 *OQR RAIA eT GOIN | 
m~—= MOON AMACKHCOCSC *COCO *O el #O ef & CNOOUWDONNODNGOOO0 
Y=H—H HEWOCTCOMO *MO *O *OCOO OS CMH eHOMINMAIMA esi QNDfatTHeHCms 


C) em OX MOUV © © © e eM & OH OLN OM 6 OM ERINMMOMNGCHMOTHOE AINA HW PS POUN EA | 
em COON Ca AOCmMO IMO 1 OMS | OMOS [ANI HI IMP OPoO;yMeHnntomnnny 
QUANG a Cm ee tm en me 


POM FAN MOA EF ON HD FM NN AOL FOAM A MIMDOOLAAH AM e B&O wm we we eh MO 
—Detret Jy meer eee ®eee eee e npeenaeee o meen ee & BPOOCOCO He atta 


LN ee ee et NN OY OY SST STL LO UN LUN DO 00 Re DO OD ON OF OR rd ad eed aad aad nd dnd end em | 


NOVO CO DOOOOO OO GOOODOOO OO OOOOOROOGODOOOOGO00000) 


Om We 
oc WW 
mio OO 


: 


TASH ODOSNMO SING DFORNO SADE ORO MINMTINOE OT OMNM GIAO Ono 
PRAAEMMRROOCOCO OC CO OG O wt et tet at at St HNN NIN NIN IOI APO A MM MM OED OS 
DO DOC Cotes ed bead tt aad ead aed at dre eed om em em ed ml = ed pnd em tid ed om od om pend ed ml md ed em pend eed rand ead mad wnt ema ed 


I | | i 
I 


=< - 


127 


COSGVWOOVEEoSGOR000 000909 DISIWIIOIS CI0ce00009090090¢C0 
NMNMNNNNNNNNNNYNNYLNNYN YNNNNWNNNNUNNMN NNMNNNNNMNMNNNNNMN 


sae dd ced ceed ced ced ed cw) rd cae) ce! ceed ed ee! ed eel ed ee rd ceed eed ted end ed ed ced ed ed eae! SP ae 


(IeX€T),F201),VY01) -DELYC I) »W(1) »f=1,.™) 


[med - 

aa a 

-_ Led eo 

(ox w~NeX = 

% > i iT] 

= ahs ont LL ome Cd 
> _~ wee o _~ 
~ oe =~ — ~ _ 
ie 7 alos = a 
— we — _ Sa od | 
q> 1 oa] bh — > WO i) 
Cc. i ml) wr LU) Y 7 
(Sy en > + Mw o o 
e za = — ae NV Mile = ont ~ 
= pee ay eas ea — “+o 7 
i N #eaICHt a st mew oo Nea o 
N YO -?>O~ 5 ROE oN =m 
= = ioe J | pal 2 meer eee (9 OO o “eo wa 
q) am ot ost — CO CO mt tet HO <a 
= c— y~moOT— “JO wet Ow © ~~ ONAL ee 
(Con Om sy eOe were Cy > oo © COG 
° ~oee m Z+N OCewrenPrenm = Ce = 
(IESE om ee Oot Cane FO OH OOon OM NOT 
rome See Stee HORDO +DOSO Weve 
vies CMMI wOd Cem KX MORRO HONDO —N Otn 


Tie ional Leo MO. eam owe eee ecm ee cue CQ) e 
aS CHEK MM Ter HHDil eOU—-$NND | OOOO HOO000D I lc 


ieee > PMO MN Ze O HMO W UW ZO ewww www ZO Fw | Se ee 
OO HC Hm mm CON Oe SW We Zs edad 
aS TN OO TO th Ot et et em ee be be be mt vO DU EEE 
mm LO mt me Yl re me mere et ee et ee ee es tt et rg ee pet Ped Pt te eee Ne pO OO 
~ONA=H=ULOXAITULSCOKMOC~OOOUEZH ECU AMYOUULOXS2WOLTOOO 
COM> hOM> MDM DOOOMCOMNMNOm EF FREMBEBBOmKOBYTUOUUWL 


N iN Meet SOW! om NN Wh ~ Oo ODO HAOWNO 
wn Co CO COCO COM An & DM QD mete 
O = ee EEE P= Pott met aA OD et NNNS 


NO TINO OROMNMA TINGE DROKMNMTINOR DOROSAMTEOFOTO Am 
TITTIES TMD MNANNUVODDDODDODOOR RR EREAERREO OOO 
ld met gemed comd peed reed cred gee rend reed ceed weed reed oe met ree) creed ed ee el ced ce) ctl ed ce to ed cet med cod cmd md coed ed =) aml pnpel) pam 


MOOCODOOCO DO ODOCODOQOSCODOOOCOCOOOOOCOOOOC0O00 sacs ad 


128 


-_ ™! oz + 
EE =—-— OF WN 
oo = ot (©) 
wr —_ Wt) +> 
tt Fee — 
Ww Tmetw 
a roe £# 
— SM oli Ql mt 
" ow>a Oo 
ot ot LL | 
xc =mwON > 
AC mG bf 
~~ while w= 
YO FeNe CN 
— osu YN 
AeA H OM 
OMnON Dee &H O© 
~~ eINNe F 
QMOr- rT) Lu | 
=a AT ee mm LW 
Dw Av o) 
CO Nt thee oO! 
eMeaa OO UW. 
aAr=—Or Uw <«f 
orm ON OW | 
Ow OK C WO 
a DYueo2> OO & 
Wu Ww uJ 
NOx<xWwOm FT + 
MDenNOerO WW ><) 
YrHtotl— + 7 
Com WU NON 
win iin WM om 
ae tI >w w+! 
Dwi I et 
OODLE - © 
em tO = 
mux<eAzur © > 
oe NS NN 
Arp eS WL 
an OC | Se 
Ui esmQ e@ ec 
AN —~X Aa YY 
omtrmtuer J Oo, = 
Wax<ao w = 
-—Orz~ Ces 
reyOuw Q -— 
No oa ZU 
- KO US ome 
NIKE i CMW 
tot (0) oot free aT Wot we 
Qe Ce wWOw 
sd ohhh wu: 
° om} 
DO on OW = 
=—Oro 5 ae & 
Iwror -—MoO-— 
AITO NSNZt 
Ce bt 
Otuod 
CTR 
Cxewe 
wic+O 
UL NLL 
od e— — 
a NG am Oo 
ac oO 
=, ae Oo 
wn 
fe of) 
ond 
©: 


0184 


eECcoe°o”g 


LS 
LS 
LS 
LS 
LS 
LS 


H 


LPE15.796 


P[2,2H)=1LPE15.7,6H 


PI2,2H) 


0186 


ie @) 
ad cmd 


0187 


GOVGOGGG00000 
IMMMNANANMNMNMNMWN) 


Se A 
| 
| 
[ | 
i 
: 
a 
imo | 
| 
s | 
x< 
— : 
@.: 
oor 
~~ 
ju ! 
= | 
od) 
C , 
Qa — 
ie 
z,.m : 
OW~ 
= oO AY 4 
He OD ~ x ; 
(or Pt “a ms 
< «KO | ww SY) - 
=) | = ~~ On | 


¥ =x | e {mu 
ad ad LL ee SO HE Oe DIO 
Iteqgil Wx Jhb 
LW Lems jf oe Ce i Ou lw 
CLOMYO—MA 


Oo OOO 
a NOT 


MMO SIN OR ODO 
poe ooooce 
ODOOODOO0O0O0 
COCCSOMOCO0CO 


a | 


129 


TOsRC,»WLBM,WCHAN,WPERO) 


5)**2-RC(1)**2) 
SQRT€32.174/(53.35*TTO)) ) 


~~ © Uaen(jmuw 


OO ee 8MOW. 1 WIN 
LY ¢OF eter enum eo 
tu O eM OMPaiFOe ce oe | 
ZZZAHMINAH ON HN HH 
poet J) et eee od re (LS) Lj om em em | 
l= NONNMT BANGIN 
CZOr OKRA H WH TOOOOOWr 
YweteTaAtasdevrvrvyrrD 
CcqzTESEurk teOWOTwWwwww-o 
PmOGOUOTeYWVACaAAaAwWS 
NAF Partd ESTP BETETEOW 


| 


IDOOCOCO QO OO eet tte tot nd | 
jooooo elolololololelejlololeo) 
ODODSCOQ@COOCO009O008O 


| | | 


ee 


130 


Cs -~ ~ 
* (‘Ss an 
ec om om 
Ee ~ = 
(an) «x hess 
- © om tl 
Ma CW <i 
UU. 
SS) 
<I 
i 
bom LL bh gq ped ~~ ~ = 
AY Daw - = ia: 
xe Ox = a ~ 
ae am eifl) o—) a 5 — 
mF Celle <1 < <a 
Nt mhNIO _ hk LL 
eLl ae el UL uw { a) 
Sw aA N NJ << 
aq <{ke oe | i = 
>e > eM “= ~—m_ + 
aN em & On Var? a) — 
wt om OE = — a 
>UNOR~C au < 
ee? eww) be bk oJ — 
AN e~-E FF — iL’ <f UJ — <x ~ 
INH add NIC NIG ~ \ _ 
EeOrI>a +e ~ tt — + _ =~ st 
ee 9A WwW om tm = z ~~ 
TAWA CON ao =—m x ~ x 
adnjCemre ay) AW , | < | 
e<t Oem Wa ny : = —_ = 
CO tay KB ZO faa | ac N _ uN 
hm OO wet la \ _ ~ ~ 
O<q © esi~ rw) wn bo =~ bee 
Ode DIY NY NEY — mt — ~~ 
Pm e}O>< e« om Of ~~ ~“ = ~ > 
fj od Li ed Som —_ se = me ~ = 
ad Pe tS) Ne N= -_ <I _ > 
eM e— OC a ml Pt a) co) LL wt ~ 
“ elf emi ew HAY mu) — ad — NS 
om) mw we Y) Soc Mac = < Caml + 
aay OxeMUo aa fog < ! < + 
<< e0O Ud wt Ni (_~ tL - Li - 
LL OMe i NS Tad) — tO a Zz an ! 
AJuLY xKOOU oom eee pnt wes @ < ~— Im ey 
aK eo & en +O om de | t om {—-— — 
wl eee em Oo OH “ow A Ome Om <a ww om S 
QAOOO= a“ Cw+ : wt NO ON MUL UL) ts ef Ore 
CK eNotes O ot > om - Send on Ow mL wwe Jl ww LO at, 
CNS www [oar eRe eM Khe w | er ; Of etl. x OIL ja oe 
RY eAtNadun OO) Cherie Mew MY WY ~~ ow olan t. od eo 
qI eyvy<at eucd Oot ee te ee ae | Ou ! ~~ UL <I UL see Ow 
ee ele #~nN Mad aq ace fe Mus matt = aoe Ad om 23 
MeN JOdw TAWA FU RU ed MNS azcewe Neen me WI 
lL MA dteaoee MH eMlLy Wo Wl evn tf ~ eos we KIO et © ow a 
We So wee SHON NE NHI | Ol Ane TRRHO HOt Mt ti 
Zeer S2>CCc!l e— NTH HORNON ARN 8 Ne NHR SNK es Py TOMI m=O 
tl eet el) NO | KOS etO en eee eet ed (Deer HE POO LUI mms OD ie 
bh <I et om re Ce eT et ee OO eee el ae Mme KENNEL HD = OM | 
MY MNS MN ZCC Hm wee NN {<< WV I —K fm BKB Om AY} 
ee NKZAOMAG © 6 eel MNANCV AON HO Oma OO | O-—OO wi yo dagen INO O-) 
MM euler YR ZTOCN SM < che de be i IM we | mr L emt tO 27S Coen | 
tre eSNAWS I IH wh LL ee he Keeler Le MiP ep 2 ae we 


SS Hoodies ein no imei aolaxonneceeoCEubencojo! 
LPOSYN © MOOLUO=N J ON TONTINOWEZHOO=O~ OOF ANNOORHONNGO: 


=AANCANO ei 
© in NO = Oa mn § Wt OF © O 
; om © o- © Oo O- ast '§ OOF Oe O et! 
aa roa) om mM moO mM MO | MAM MMA MM mM; 


ATFINOE DAOANMHTOOMOPROMNATNOM BDHOANMTFNOM ODORrN 
OOF OVO Ot tt eh ht ot eH MONONA NOAA M MMMM MMS Ti 
See eee aoe OOOQOCOOO0OO OCOQOQVDOO9O000000 
elelelelelelelelelelealelelelnjleleleleulelelelelele) OODQVIOOO9NO0C0900O 


00Cc1 
00C2 


i} 


' 
| 
1 ! 
i os ' Sat eas he 0 


131 


a % 
Im = 
=m ss <I 
= am SL) tmnt -_ 
lent aa) Une = 
ae tee +> <I - ~ 
KK AL — ote 
1 % ~ 
~~ % <f eo} °° ond 
<< =O) ad al WY "awl ' 
ae mets SN OOS ran) < | 
~N a+} LA WN = - , 
~~ wee Lee eH BH tet | 
NT aa) - aN Ot “—< 
mY oO ~~ Wi Wi Le oe 5 
—i! NM LM ++ +4 =) Om 
=m CO mf ee ee <I LL om 
mat % te +o - n~ 
>< wee wee -~— += Am NC tome oN) | 
1Lmwn =~ ~~ ONS t= 
pes CU me wer oss on | 
A+ + + N+ *H H+ z= Qe ‘ 
Sao % % HH HH oO Ue | 
<— mom om tC) eeee st Ost med Ginx! 
= OS = =~ tts = z2 bn Vere 
we eee — tat ~~ AQ Li WJ Lys eT eo) i 
=~ >< > ~nN anne ++ ++ x aM 
me | (-_ te > ata etned ate atme e°e ee Lu ie : 
oN AN NN #4 ANOS 00 OC = — ' 
* mime SU SP O Ne» er ee ee ee ~~ ONS _ S$ =m ~ i 
+ mo) wee mee wer ~~ WO x< <>< MON MAR ou ml i 
= x | xX “Nw oe ee | “+ + Na ~ ow I 
pete ae KT ome ee N™N NO tn = OS | |=) uw | tomtom pee ON 
ete Yom ee —~ He MK KK Wil Wi a> oer NH 
(To - — mt HHH H ee ee oO om wt «OO 
— =—~AHNNYH™ ww PC) Jem ee ee NN NAM = No eed LL. oat best 
ag Oe eee ee ~N NN ce qo 
ewe em WK mM ENR COST IN NN NNMONN Ww +o ewe IN EN 
‘Ci ewww Ce ae eT eel ialial) HH HBHMAMM mH ee we Om; 
Ol NNN te SH MeN mee HON MEN eee. o CN me ee Ot Re 
Ne LO~n~—t Am AUK HENNNN HH VSFZ Z2ZMwrT Hw ° +H HROTR NO 
NAM MY me eee 2 TOW ee x Nt er et 
‘CONE EHINM OOD Ne WHR MRK DD LHe + EWNONMONOT NUD Whew fae 
SENSI we OO OF OW AINMGOST WM DAN MT ee oP e eW we tt od (>) fs 
mom YI 8 Od | ONE Nwww we t NZX ZAOTANMMN BPO at etwe Sp ww 
=—TOwOOO~Kll OM Cl #eULMNNNNNN Fa a eK oe fe EC Ike HO 
KHAO HN UN tee EY ETON SH SE EHMMO Ow SOW OP >ddO iO tl 
~~) So mm mm te MSY HNN DDETZ ee ee ol mmm ww -—-D> Wto>reaere 
HOI NM FOO I I mH Ree SH BPN DDN AION Cm ZN ee Oe 
Rm eet Cm cme HM PE ee ttt om I] ee I we 
Ce CK DK DEON eet he ete Ue IOUS HM I mm em SD FeO Th ee eet eee 


WWM OODOO wwe wrt NESS Nh SAW PD we we wer et dee et et 
ADKNNNNMNOUNONEA S$ DDIDDIZZZZwewwww hb Utero COecOOd*D awww YX 
MFANOOHDOCOHNONWLONNMNWN MN WWW WL > eee ELL COOL“ ROO OOK HOO 


NS 
oat eed 


CAC 


MOTOOM TOOANMTUN 
oo oo oe AL SAL AL aL S 
COOOVD9Q9O0000000 
elelelelelelejelelelelele) 


essa 


at 
=! 0 co MONS O 
N ond om LAINCUNS 
(oa) AC AMAA G 
OM DRAOAMHNMHSTNAOP DHOHINMOITW 
WINUVIN ODDO ODDO OOP RR her 
O0O000°0 OOodo000O0O 
OOOVOOVOODO DOGO O0O00O0O0 
| 


132 


Se ee 


aon 
fa Taalaa) 
mene 


oo oa 
al Bn 
oo 
oo 


0076 


iS 

aa) 

faa) 

i 
anal ae 
DDODODDDD 
OQO009000 
OVOO00QO0°0 


— ~ 
aot wy 
uw uu 
= a 
< <a 
a e 
* Bg 
fa Cs Ca mn 
+ + WY 
% + @& fod 
_ — —_ uN 
-~ Y CA Oe ee 
Lal <I a ™“ = 
— 3 & = ie 
m= <C lo NI Ce ca 
me TN <I as G 
be OS OO On od 
Yr eM > + a | ] 
AtTO@Q = wet FO eet = 
O #mt ~~ ~~ = 
“a = of WY “ACa fad 
mo OO e = = iw 
i emo <I aIm+ + 
? Oo as — & ~ t~ em Ee 
WN Tae OMe HK MH NAM HM MS 
no OM 1 OY tn Gn <2. Ot a 
: ee AES moO 2 Gey a <a 
OU © oct +  ¢eu U WM WN 
CO ome | | a? LL Nm =) =o oe 
Oeimmnm iy Ct wT AMA a < - 
ll Oecd NINH tt 
Cs | | © Sd ed © © 


tN LNRM ae ee DD 


ALON)-ALFALCT)V/OXONI-XOT) D+ CALFALCTY—ALFAI(M) )/ 


~ _ 
~~ ae 
—_ st 
—— —_ 
«< x< 
| \ 
-_ -= 
C. uw 
om — 
x< x< 
—_ — ~ 
~ ~ = 
~ = — 
ae ~ « 
= Si ~N 
os a) AN 
— — $e 
<I < + 
a uu -~ 
ond i= = 
<I <= md 
{ {= ~~ 
nm mm eee eS ml 


T il — il = Il he H 3% © 
mA Nm SiN zHat ey + 
HAN Mee tNIN eZNOAWIO 
mee Nome MEI} | DM 
1x We} MIKNAS 2 


COW ef ele Oe CH OR PST AC ee Nett OO 1 OPO UO Oe 
(DODO wee wee wee Tb bt KT et be ee Of me LL ee be LL tO wer be b> 
im WN [0 ee © Gn > Caer Gogh) © © Cn © ca ee Dl Oe A ol 
OUCAWUL WO WOwToUUCHOWCC HW WUKOU I KO Mee GOW 
Oe Odea qodeodOodeqOodyodoOtZ=OO™=O~OO KF AMNOOMW 


Ye) CTO NN FTF ODD O &N 
ina) Mr FTF FT TEMS HW TIN 
aa) Am MN MH MAM MM MM 


! 


=! 
Ove COW +r 0 
r rn Law wv UN 
O Mme fap Tae) Oo MO 


IOROSANMTNONOTOANMTMAORDAOGNM EIN OR DRORNM TIN 
DOAMGTAFARAPHRMOOOSC OG OOOO Met eteteadteat tod tet WON 
DDD DOO OOD OO O ew wr ed at rat ed red ead cant mtd ead et eed pend onl ed emed pamd rand oamed ped eed goed ld 


CODDODOOCOCOOCDOCOCOCCOCOCCO0CCOCOCO ODOOBCOD00O 


| 


133 


C4 
(e4 
mw 5 
e 
io 4 
<q 
tf 
WY -~ * % 
fed © * ~ °° 
& | ite ll ~~ r -_ 
”Y Sd (sa) r — t 
Ni Ww od al % — e 
o b= fod + *% 1 — © 
4) ad t = - * os UL 
<a o te q — -~ -_ Lal 
eS c= uw © — w Lf T] 
fond = mam ~ <q fo a) »< {en 
— — + > af ad 1 vT 
~ Us fon) <{ Ge Om — aa 
Ww b= + (aby) | * ww <b 
_— ~ ea e = ~ > 4 Tx 
~ = te NSN <I = — ~ 4 
Ce- © = — OW oS: ~ + - 
<I fod be UL iL ~~ (ag a ox 
“Om ~ ee tr N - + ot 
wa = + an Ww ~~ t 4 ~ 4 
— e a om, I >< o< * ~ <q oo 
LL —_—_ -— c Cm > +> + fo @) % Ss 
e272 Oe + adam Wtf m ~ eo [a4 whe © © 6 
—Y AO aN] 1 oe Li. — N La th UNA 
Se wt * = — eH wet we + — aA 
rc <w Po << | fa) Om wD + ~ Team = 
x le -~ Ort Oe t WwW aod -~ OO . Qmannan 
~s —e aa wt FT m= <x Qaee Loom a eons ood (\] CC) 
anp O-z —_ ~O oq ££ +> ke of — en) DULY meee eet sas 
It HO fod tw NK a TI WO ow led Ih x x & 
i eet) % <~ ON wu NOt * tet MOF] I I 
ule Q~ a Oo Tm NY , ef mm a ~~ Ct eon 
NW fe en ~~ we et NN ¢ ee roa) WW ems 
Keo We Cc %t % © | = >< O< 1oO!l co fo 4 ~ OCI OO 
It NO + -t+ TSE | ++ SCTE + HO tO) wee eet er 
we ef) ais as eet Zt MAI ae -_ Ost Tu HH 
A1Yt mw aa) om = Wi hs m ASHP O Ree 
~~) Om ome —_— e mw (49 % A et wee oe aw ST en & eQen = 
=z we Slee a ee a) OM RF wr Fm ae & OC ont SC tt Fea es ee 
Oe KO —_ coo ems <t & Ad + WNIT EE Pm AQTO GS Nw 
—a de NN aeestestA+tOd Ge H I NREO-vod NOAH em HEEL 
u> ww wen tHe eO~OINOm With «A +0 me FS COO M+ + + 
(> fel{e? DM COUOU Hw ef YOY — > eNMT YL ON Mat © CONN tT em oem 
Wap wa POR © eNO he ~—>PEr iene NUNN + FPA IA 


Z2aq Zaz ~e eee eT efORMCTQUMERZZDIAO~ (OC PK tate em ON ewww 
<_ 


DPBONRKF® M4 NATrnnnnnrwx Oe | WWW IN OFT NN 1) a eet ee 
CD e%# ZOCH emt 2 8 8 ee OrNA 9O MK K> HU NaeTanaOoan nls we rwcemn wu 
CSIWAES YN tsetse tweet tt WES I OO eh Ee i EHO 
CeQEeTSBOOMH Ht Wow woe Wo WNOS eH DW Hm ee ey wr FET 
De We Ogi leon NOY TOUWOW WWW Z ON CY~OLe rei eOmDID 
ACYOSPOQOOMAWUILLILU OA Om >< <M EET OK. DOE KL ah et eT oes ELL EMI 
ome Or f ' tl 
o) a 0 DO lawn ooo 
io) oO \ (@) Ora | oma mel gxad omad omad {NJ 
t+ wy + wer regress 
i 
| | 
. 
=~ NM SOOKE DKROHNMTNOEDOOAMNMOFSNOF DMO NMSA GEODROANMSKN 
DF O00 QO0C0OQ0 Ada tet ett MRNANNAINNANNGQMM MMMMaAMMMHS SEs 
= CO BDQO0COCOCOCCCSCOCOCSOOSOOGoeoeooe0goo ecoqnpegoocoooqo0o9coe 


oo ses ae eCooooo0o0goo copeQq0oqo°;“eoodae 


' ' H : 


ee Pe ooo eines a L = ee 





uw 
— 
<I 
fod 
x 
Cc 
— —_ 
CN Ww 
e 
(oy (AG 
e- 
Lt 
wr Ww ~~ 
e e X 
oO cmd rd 
va i ew 
Lu tw 
UN WY (ob 
~ WU ~ =z 
eae ok N o 
Ww WO need 
> 2 oy 
mee Vb mame 
— x< = x 
wn << ow 
e - YY _ 
ot QO «a Vals a} 
«> iw = zw 
LNs YY G Ow 
NSwn = wow LU 
ee D> | a 
° Neto GW 4 Ux 
A) Weert I WwW qe 
~ ({>u aw 2 fo Ae 
~ HN ! az wey 
-_ >a e -QH vB —_ 
+ = Omemuw =< 
a VY maZEtwb Oe 
=< + z WHat DaD i= 
! on) ~~ eX >Q ww wea 
~s tf co b AOD o EB oT te nw 
ns ~ TAL oN Dm WN Gond wae 
—~> A) Tc CO Ge Se = mute a Ww 
xD = ~~ Mm eo CO zD eZDTA =z mm O 
~—~+ WY wm + NC © NY Mequer | BW o=E — he < 
HOOK We eo MO © —+ Ae est eo EZO . Te 
KENNET HEN fe ® CIN loz zOcd—¥ {om 
-=DfND wfwDM “nN NN 2 OTHWIDITwe ww 
NWN aNM NWS Ou uw DoS ADDTMNWLWY | aw 
—++O8 EF +B mnt Ww MWwWN Lids x TS > I 4 | 
BOING | aN Tee S44\-« “” pee Yio | 
£27900 auc ao OC Nt ~ Wear 
+Deww —WO O+ 1 wEE ONTYTOwe «OT® | 
—NONL NCOs be = TDD NNONMYS weraroro 
STe*NSE SE eM LON MN eDAN of Png To Ow srr fo 


wet Ter OC w f we NE NS ON we we rm ot OX OiN OM Oo ON OLL 

SHE—MMOMDMVUT et para tu tt teow Ow OMOnyo~vo 

—Yf Bq Bae gpru Dit wewesgng —- &— KK KRY F&F FZ. 
iPrTrwron tO ONRQa~~~— ~~ Z2Uawawawdanqnuawdg: 
it ESOL ILS QW EOE ONE eR Te Te EEE Se TD 
SE I LL OF es OF OO OY OR OO 
DILL OManKbMOOUaAOQAAaAAAAULAMOWOMOeKOSNCOCOWZ 
IN Em FOOBZOs@2OUOH>S>Y)> FEET BU ZBLEL EL | Fu sucw 


= 
Oo ao oN ro ON N t+ O fi OO WO 
N N mm AM MS. qs N t+ ts) + WwW 
t+ tr +t ter =e tT + esr - + 


| i 
; § 

OMOROANMSNOEDOO QAO INOr JRO mC EIA lon omnoe 
TES TOMUIUD AMINO OD ODOOWDO GO mR RR Mm RR HPO 


DOOWTDDAGOOC SSCS OGOCAVI COSC NHOGoeeo0o0G©o |\CcOo0900D 
COD 9COACSC OSC Ooo AMO0COeeoeegao0o a 1099000 
| t 


1 


| 


135 


SRN SR eee Oo ere petite 


+WPFER2 gWPERLyHEg WU eDHEDX,SyDSDX1, 
) WPER2(10) pWPERL(10),HE(10), 


, N 
© 
! — 
' o 
~ Ne) 
a = 
C= = aaa 
=O so Cd 
1 Sed ) “a =~» = 
=a = f= CO a = 
AS: ~ Cc e+ font om om fatto = °5{——_ 
=~ = — cVO0O weer bed om te ten Ee ee eee) ; 
: oe oa ~ —™ ame DK wae font SK eee oe ~ ome > eee Bt Sl wee om 
|>< xS> tty i~— am |< (—— oe ~~ | xe | <a 
ie Oe — ~~ om pot comm ome oe 1 oer) at, =m f ppm | 
plan) So} = mm ae Zee eK KT OOK 
ee a= (aN) ef AJ | mK mee Mw OK we | sae | me mee MK me | 
<X<E as ea SR mH KHZ KE KS HE KL HZ 
Lo = uu Cet wrt PRK UW HN UH PNW 
OU mK a. WO eee er el 
lea ~C x< Ow EF KR Oe ¥YOMOF FOKMKUO RH UO © 
‘Lie m @ = UL, et PR me OE pet CO wee ese ete et me et (1) ee om et ee & Go! _~ 
lan Om — bem QL weer me eee LL SO ee et tee FA met ee LL tt ee et C) = te) 
I vend ed Co) X ~—ZANeHWK Te Qne Kx NAD Kew >< — LW 
Se Lm Sp (WU) Sp CRW WOW+=4OOM W eMWOWw+ KON > ez Wat 
[aA an = 2 o WiemOo oD ten |{ MO B@CLlWITt-~ONoOaA Ox oN Y) 
WO 0 =~ ~ EN Ow FO!) SON SOF FA /TOKNND+EMN NH te=™N 
2uaoznrz om | | ENKHORKV Il — | MN MNP KO | ONG EO ES Fe 
'mee OOOO eee rT} et wee Oe LU me tt ee ew mt Zt ee 
jILCMNwS “Leno = © COO met LL) Ome ff mee Ome YY mt | LU me |] ee ee ) II IDo72NYe 
joo zo WN Ome RE | LN KON 9 | Ow FEIN RMN NS | TMOMIZ NIN I 
1X -4IWOTCWUWMN + K+ | ewe KES I Om If me KS 1 OO I Kt Se I 
Buedtawerni>o —f60 —-— VW meer 1) gd LL mo fH Qe " | | ee fe  -  o P o) 


PaWeLCOr~aOgtlsnOn VUUSZLwrwwNKOSZU www WMO I} L~OO~u 
NALOOVWZNVOR®OZO YPM! KHNDO MN DOWNY *XHODBMVOOVNZOZ DOOKU 


Qa ae 


S Nt O co oO N © DO N 
© oo oO Oo exe — wt RN 
~ Sf ~~ ~- ~~» ee 


aA NO TWOP BOANM TAO OROANM SNORE DDOMNM TF NOMOROMNG TU Oi 
iO OG COOCCOCO Mette at Hat ANNAN NNNNNNMOMOMMOMOMMOMOS, oT TS 
©O CO DOCDDOVCDDOOCODOCOCOC DC OCD DODCOQDOCOOOCCOOCOCOOOCOOOO 
Oo OD BDOCOVOOCOQVOSCOCOCOOSSOSVS OONQMOSOOMOOCOMVOOOOOOO C9090 


| ‘4 i |b tla lel aa Baul oe a 








136 


~~ 

oo 

x 

SCS 

WY 

OQ 
wv 
* 
% = 
we Ww 
fad 2 ~ 
fa ! oom 
+ Y — 
~’ = 
wn >< 
— (ea) 
is @) ST 
+ ‘am 
foe) x< be 
% eS = 
% tL! — 
LL x oe 
o& (aa) a 
Yy -_ o 
% Lf om 
= oi 
t Df ww 
— = we 
@ Ww aa) 
+ az Ow 
N Ww oT 
+ rr AO 
& Uw o 
ih Rat 
ce SS ate] 
log = Cue 
% mw Ow 
mm x eo 
oe Nm 

i Ono @ buss 
H t+h—N WW Ow 

wen 2S ih 
Woe wet fb 
ANA J tem 
een WG r 


j NN eS TIN wo 

1 MAO CUM & ety 

‘KN RNR peer 

Hm Sm oN jf) mes oy 

HOO met Owe LM oc 
Omrt~O} -oOrF rz 
j-ecezewandwo Now 
TWD ome weet fre SET PS SE pe be LU 
OU ee ee Oe OO 
jwCebberCCOrOows 
UL tees FU OW FOO WwW 


| 
art Of ON 
NAN AA me 
SRN ER RR 


YK TOOHKINMATHOr OG O 
TT TOO MMNMNON O 
DOOTOQCOCOGOO0S0S 
eleleloleleleolololeleolelele) 





137 


_= 

~ e VY 
oe! —— fe @ LW P- a 
ae © % Mm OO OO OO 
x dav — om Eee ye mar are 
2) = uN — df 
w= LL WJ — ~ Ow & a. 
=e = Yo) Oo + + # 

om AO fa a] Ww ~— ~~ 
mY ed O + awn Wm WH 
be Li. Fm ae GH faa) Chi were ee * ~ ee 

Ke Oowrwe? + wm st LW = ~ 
mt CO om << & Oo memo a mal (oa) 
am we pam pe ee ~ x + + + = — 

yo al Lu WL > bene beam mm << < 
Ne RAINS ~ : % + + ee be 
>< = = oe WY Ui 4+ + + Ww Ww 

ya a rome [a 4 (oY (6) () ~] ~ 
Luce —OO >» + i) ee eee a ' 
—_ © ent ot oe > Oe no) ae 
O os md ee we) wf (2) Va OX Oe rey uN 

= DA coer — HF H+ = = 
La WO > oO No) Oo =~ -—_ = <{ <I | 
+ re om ~ mm <= Sat oS ad = 
ch ome J a A+ —_= me — a Ww Ll 

eu —~Oodt “! HN moO +s W A Nl 
WS Oaatmib= + % a} fea) [es) (og = = 
DOO awe + am H —} + + + + 
eer wre) & = os CG CJ NS CA ote = 
lO) AS Om — we bard ~~ P-a + % =" on) 
aDa at e-—O© ~ ey Om ++ j%+ # VY) —Y— 
ht ee moO ma) KA Cy) ©) & [ad a 
iS) pet LL > ened wee ( % © 2arcw ] ! 
Te Oa~wO A aa) tt (oa) Ta) 

ox <I omer HK fam HUY Om fee (cf (e¢ 74) 7) 
DTO— —eCw + © GO uy + + ac ce 
=O-+x ABO e at NOD MO Ph (MS (a9) (5.9) ~ is 
a OQ aw & Om ~~ wt Wad om on ee =~ eh 
cox<—_ Dao —~ A= Wwe OM fay Une) Se ial wd ios) 
Ou eOON —— Im +m He ~~ ££ OO OM — Te) 
an om tee we eit ee Oo + te Sd ec ao 
iq#cT- Ow~O ood er awn 27 Olt go oOo; ' 
POW aAaxa <0 =e NN Fo S22 ee HS = 
ere ww TO e ¥ -~ >A ad wom Cf [ ) — tet Ue mm OC) 
aa ole — NINN Om Wh uh We fee (4 ee Oe ~ 0 
Dae ftom ToO ° we om he Hm fe wel & Oo + & € THN WMA 
>eM SPOAN tH e oN miltte NW Wea OO mt =~ =~ = HY ao 
“—LMNO eet ewe DOO NAR K Oe SSH DINO ON NN Ww ~~ 
ee eG rt Nie fod be wee ah 2) |Z elf uy —t rr 2 fe 
wate eOHNOow eNO NOM ew DNew lk THD meso ye OF OM fF WN Ae = ow 
Cems 2 MOO eH DEFOE NS —NOO Oo 2a OD DO GCG eet OS OEY 
CONY eee YO PT Hae scnet eH mS Pt) uy © oe i SS Go 
Yee DON #HANNKSqawranee I-O=- ~ aL ae et BR Oe Rl 

HOR DP Hs ew FS 6 OO Hh ww tO pe IO mie COC ame we we UU COWL) 
WWI we TOR OAH Meee YD Wwf © OMNaIM TFT Wan nN NI | 
Ze MM ZAOOZHH WORD RHO DOM NOCH NNO eD DO ODO | HONONS | 
LULL COO Qt om OE et Lee om Ip Ob II Wwe Ol O ee il i || enon l oe Cenk oe hemo © 
fem MICO mm owe WO SS llhe>t il-_— ~rmt ONIN 0 OQ We eR tet NTE) tee 
i © FON LOZMYLOO I er=— Ieee Te OOZ | Fld oe 2. ds ae 
TOF % ZFOANOS HH WAN I ee et et mw en Hh we I] Ze | Om PUT eee ONO CYA 

qt ft Ce Ieee 


Creede EOE NW Herter Pm Zeer Z eZee te Rete Re EES 
DAO We AWC O I Ome UD WU at ee LULL ee ome ae ee UU UL CL LOW OL lL © 
MAX MO AOC MONOD ODDEDZKAlLOMN HOLL YOM FUN N14 Nt SN OONIONINIO 


CCD SW S ie a OO COCO = I at 


OaeN OAINM - gE co eo MM 
HO WNTOD O NN St OWO 
: oOo Oooo #« me et eee A) 
| | 
| | | | | 
! { } 
| 
| | | | 
+ PNT aa) + Or oronnm sino OromNm NO~ © GF OC aMsrnoro 
) olor DOO COD et Ht tet im HMHVNNUNNAINANION (NEMO MAM—MMOMAMM 
oO oo CQOCOCOQCDOCOCOQCO0O DOOGDOOCOGDOCCO |0 010 OG000000 
ro) Co | OQWOOCOQOO0OQO000 CDOQOCOCO00O !O aw CODOODOO00 
= | | — a Dia 





138 


29)9/(1X2(3)-X2(27))) 


~ 

ae 

ae 

ipa) 

— 
om LU ~- 
~~ T ~ 
o~ | _— 
Nem faa) 
~~ ww ~ ~— — 
FY = Aw N 
— | om tlw a 
Na Om IT 
ww (OV a om ww ~_ 
Ne NN Os peg 
xKWY — wm — 
| ~~ ! Nee (oN) 
on + om xT x 
Cm fF ww = 
—em er COUN ~ 
CU CI om om 3 on on 
Ss ed annl Don EE ol 98) 
am Nee Nee Oe we 
mem fmm QL WLW 
—Nt—O} i«K of ot 
CN, wee meee CV: eet ome en tee on ae 
KINI KYANICON. LOS 
emo em AJL ee ALL 
NONIAST «MIN ME 
Car ET Oe ee eT ee er el 
NNANYN + Ne 
Mw OT) Maree TO on ee 
Ne Nwmeate em wm 
met | ome SM) ew 
awn WO~ WO 
~C~K~O/}TIT Iz 
NNN NNO wan 


NI $6 UN LU LAL 


we NU\—UN\T «© ew TT ec 
V2 OY ew OON~O 
~OCOC~O tt tt tho>< tt tt 


Moth TE UE Ep eommm | moe LL) 

mer ee tee eer See EK OK DE wee OK at CY 
MxM KK KOCANCAkeD 
COOCOWWW<WLZeO 


SooooeaoaGc 


Oe 


TFOCaANMTN SC 
ayes 
olelolelelalele) 
elolelelelolole) 


AQOdus 


N 
N 
Va! 


~DORhOex 
me a al a a) 
Oo00o 
Ooo00oceo 


139 


ean ee 





CaN] 

-— o 
aq — o 
Ndwe OF a= =~ 1 fot) 
-Oae = OC eO = ¢ 
ee eh S|) wi me ~ p= 
NN e NN ew O~ a za 
Tay feOtn = 
>) >OW~- l oe \ 
ef] e et ONE im a \ 
Ait mw >T o pox = - | 
= Ca e-— ~ — % 
ele taf eO an) oO | ae: 
NE wb Ome 6 > | = 
7A. NWHoOo~ ke + -_ | z 
Sew OC Rw HY Ww ms ie tod 
AWO OC eNw wy @ Oo” — 
Ne © Mea aw — = x 
I-eu OOWk-O + NH SE 
> ef eet eocl o = IN Ww 
YY awe Ym ~ = = >t+ -—T (oe) 
NM ee OCCaAO om Ww baad —%e mS od 
— ONL atte el end —_ ae a No ew HH = 
Ode ww ewe w ke —_— —m WN mam wa = 
NI) atom Oct i! Uy ome x ~~ TH ~ 
Cee KOUY tk [a @ CC ate } aN 86-4 uw 
CNC COmtNAewW (am) QO ww ~ wf Ue e 
ayo Yew ONO + tet NN z NS ee mo wn 
KYNDM QOD e = ~ «MM —_ Ie Ne Uf 
Mee e «OD em ag) ry 1) N en Nw Ww 2 
aoadt OoO~-Oe [ed ce NW — m+ + > Lad 
CS Stet ee ! ! ed ~ — he et uN 
> ar a 1 lea > 4 = mW AN a Pe ten | 
Comet KAY ENO a oc <x ~ —— mt a WY 
WYO Orxek—-E fad Cg aw = — Yet tee 
ZT ew WwW! eOD ew — ~ YN” — ee a 7 | (ea) (&) 
Cre Tm HO ~N Nn em (oN) ~~ mA Aid Zz «a 
eee OOOO -~ -~ < N -—2ft a> qa 
Wa mn eat alt oo ce QO wy — >Sent fe f— wee Ele 4 (5 
Tet ewww iad gC aw Ww ~ ell WN uved 
em 8 ONY wWO a ce NN me) Nano an, ou 
meat aD ef ! | << ] Ma Tee HCC ts x 
IXxC w emQw~ om =~ ee ~ mr We nznwn 
reed WeO ow Ph mee eo = z —- -=—WU0 On-O me =O 
uN e TOcaiesT Ain ww MOM = ~ = «KOO Oe = 
Pa E emwod ad wer 1) aN — MONE NeW Name 

e— OA eee N 6) ~+ + NW bee Nie NOC et mt Alu <f 
NeW SDrwo = Am om mee Um dued+O~n we emt LW 
xO w> eke 6 an Mm AM AN OF RK Kt D CEN NSN 2 
TO emt! em pew co o~- uw ew Ad we WOWZemHe KW DY = 
me & Wc mO oS ~ wN EN NSN NN ee oN Creme # OD NK wl 
CeWSE-Oaan da = 3c NT a «a i Li RM wre Yow He ewe om (mY) 
KEIO ella OrF- CO * ema FE F&F COL Ya | ZOZ+tDeowONer att WC 

bem CUD OO wee ON LL et Tr m—~eM™Miy W W ww eo eNOS O#+ Dew Os ONIN 
WU ey Ae~—m~G ~- Ms Om aod WAN C= HNO HANNO | NaIMHO 
ZNNXKZELY eS “as SUNN mel mw tt at Ot oO} Kew mM ON 6 
=m eH eC © eYim ome ee New | em em rem Wee mY er OQ INH OPOMWaAoONON 
He ZK mem CO ON ST <a et A or Ce Se ee we mH ON EEN I HW) om =) fee ee 
DK YLNOODKHHBOZO je e>vye~e ww ~ZkKxKO mK wee J be) ome |! MEZZO be 
Ce OY Set HtO SHO fh SHEN ER HNIC NON MEO DBA ONO IE 1 met be ee oe et FEL 
YD CU ww et we DUNN 8 ONY Thi he TR REN | oo em mt et OO we we 0 wee wt TE 
CS ONAN S NSU YTEOCNT Hw Re PU ZO ZN OO 
DNE YR Dea NOOZZ i AHOUWOWOWOOMO I IMZOd gear wOdeLOurvOad 
AX eVLOEY MAO Os OOOCO—=DUOUANYDVOAOOZZO~OrANLnYNY& dy wOm EU BW 

aA OO ANAT md 

Oo = MS fe] © | Omi N 
~~ ee oe) ro) a CON SN 
=I N ATNOM SROHINOIEOHOPDKROAN ATMO ODROMINMTNOM OD: 
(S) (a See oosaea aot ote YANN NAAN ANNUM MOMOHMOMM MOM 
ps (o) CDOODOOCOOCOOOWOOCOOC0OOD DOOD OOOCOCOQOCOCOO0OO 
(>) 


c COOSOCOOOSOSO99000000 OOO000 OVCoOQgo000000 


{ 


ae 1 ee | co, i 


140 


HHL 
- eH 
om am 0” 
om om He 
~ [Kiev i 
=~ wane weer (_) ee ee 
~ CAINS Ooo CoO 
aa) ><I e AN NAN 
om | ~ bet wee H + ' 
Nas roa) NH HH 4% 
xo << fe > om WwW Wi 
ia) ! -~ LL) CA ome ++ ++ | 
~ ~ i CH we o°@ ee ' 
co + wt ~ ~~ HAI tr ++ 
—— — Ww NJ NN NO 
N= N ale C-k ~N ONAN 
x NN x % Wm iu te te 
XN ~ WO tN~ HH tH 
mom em + —JINY HO Sin 
one | ~~ _ tot fe oeee zgva Pas Pe 
Am +m ~ wi twr9res Www wu 
were (0 er ee Leal Wm | NN ++ + + 
(aie a Ala [ot ee o Qa Tem ee te ee ee | 
KN KY + # GS eon — om a oe OO OO 
N 1 1 x< tot + + x mA tea t+ ~nN ORNS 
Sy Nan Wt ~ =~ WY) H+ ~— + ANNAN eH RH 
Cl on KUN KY we OC RH Re eN a oN cote me Cie tee 
— were eee Oe ee ee ~~ NOS Wit LULL) 
oO N+ OO Ne KN ON we Ne CO) wee tae wet tae + + + + 
alt eer ao ome it om KN HE me IAI e° ee 
Oo - mame | ON eet et et OK OK oC KOK NN NO 
NON UKHO —N ~ONMQ NX BNO Hm HHH NN wn i 
-~ | @& OCK—- X&€K Oo es OF $0 ¢ wr rnne HH Ht 
tons e 8 ¢t) ACAe NO BHONO wo FN Nwwmnne st wh Ht 
~~ ad - ~N1LOW IHW OF OHNO eRe BNO He an <tin 
> ~ | mo ON we re Nee] et ret CO I OO eee Wi LL) 
omen amt tm + momma Owe Kedar tm + OOO OD TS ++ ++ 
OV Heo om AMLUINGI GA TV ere eee MAIN OS OH RIM + + + oY > om la 
aD ee eer LL OO wr (if) © ewe YY) Owe eM cl NOD WJ emma >} DWiw seer LZ 
e eee Om Ow CCN RKO eHtNSAINHOMH ONwe OTe ANON MeaWw wu 
see ii OCOU HX a Oe MN ee SN tee hw EN +t +H HOH + 
Wher ypuss DS mm |] Mr WHNNNYL MOOD JNO eS DTM HH He 8 ew eM 
eK 0 ER IANO RNS HN AST NR NALA LN DE BH eae I 
~ Ne terete RR RR RN RR RR SSS SS HDD HH 
SOOM HAT NORNNS ANOS met ttt) et OD ete err HH HE (NNO eR mmm 
> ~WDWi CF ~-QO00OKOQ BINNMe fe SNe te STS Sener wn mae tae See eer eee () 


TOTZODLZOLNNNS NNO OL rh Ht Oe eee We ODDS SSLOOAKK YX AXO 
~COdtA die ONOOC OOONNCY=$& a laa a ili 


=! 
— 


oO 
aN] 


CANAL OF OHOAN 
vires Tr nnn 


ral 


OCOOTOOCOCOOCCOS 
elelelolelojselolelololela) 





= 


Cc = fas} 
roa) ron) oa] 
AFTNOEOAROs BMY 
Cee ee 000 
ooo = 5 Coo 
ooo 


141 


oO © 
~s uy 


DOS DROSINMSINORONROAN 


OOD DOREEERARRAR OOO 


Ooocoo 
OQ0000 


COOOOCOO0CSOS 
COOOQOCOOO000 


2CT)sRIACTIRIZCLT) ,2I4S( 0) ,RICT),YR(T) 


) 
) 
Z2(NI-BETA2Z(T)IV/0CX2(N)-X201) 0+ (BETA2Z(0T)-BETA2 


-~ -~ 
~ un 
— — \ 
b oa 
uu US my : 
Mm QO wwe } 
oS Qa an 
+ + <x 
NS fa) eo cay 7 
+ + om er NW 
* c : e® 4 OC wwe 
a ~ { ) AUN 
~~ — — fad w rd 
_ - a eo x “eo NN 
. —) > > — —_ ~~ 
% <q aq “NM ~~ ee 
tt ian) + + ~ — mt 
-m x nN cad = Nm OO were 
i a a ~~ 6 ~— cm co ww Xe NN 
o or” mor be we) Ona xr dad 
=~ ® OO ~~ a + +~ | 1 be 
ia (SJ 0 NND (ej -~ am Com — Lu Ww ~ 
— oma un ome <auer WA al be ht = Om j —_ 
—_ OO rs S oY) pm 2h od & — wer (Af ee ca {4 i — 
Ea eG) OOo WG mL te & mu fa 4 LN CC -—- qa ae 
~~ TNOC ANN om Oa—o Oo eM = Zz Nw ww =NU beam Ne 
( = sO ewe Om war eer tot A eos & <q <q o> +> = ww LL) om, {uw 
ml eee eoCCc — & OO a od OD eWQje % — te it—_—— am ANG Os pe >< 
=A POOCCDO b= Zon = Tso HH MM AINM OO MO dre_ we LC) 
—stHNOO em ow >Ore SAR NS DOD — CAS we we ee ~~ Oe 
{ Oe tated ee em zctN (Um NBD OC 2 OK dad a Wa + wel 
CUA m ew ew el DO mM mmm ELM MIA + of +Or -P FEF OO Vote a ae a) 
HICDMOOOKXK SN COL Wetter et mes i ek Ot ik ON a Ce we em eel OC) 
‘I> CO | NOD ot ewww etYee | JW © W MYO @ DM WIN (| ome ot em 3 
Im *OKe soe AON IONNRMR MO tC ODO. TI NONONW men | ame he 
mf SZ mK I er OW I he Be {I " " (LOND LU rm Aj tmt 0 
Ch WDE MOD tee SDI Oe FE mr wee cme me Om eet tt CN et NY beet er ee wx NT 
CINKOO! I 1 me FF -—HUNCNOoereZen OM OC ww Se wZMMS Re | | 


—WWZOZMDOTVOS WHOM RZ Net FE EHOW O~ OK LNONONHOONHAONAN It I! 
mat bm fame ta ae ee a me ST be frm Pe mee at tet emt NY ce TL Ebel pe bee we <The eb eT pee | Re NO 

Dw Owe Oem New OO wt te 
Speke pester. sees len lee a WOWCOMOW Oto 
be et Zt et tt et te LL EO Eh OR DODOODOON=OPDODNYWOHOZZO~Or 2 
: ome 


Q POCCOOMOOCO om HN MM SF NOMEN OM t 
a HNDODPONONS = uw | N NON N 
—“ ado oad = ol = N { N ww AW N 


| | | | 
| 
i | ‘ 
‘ \ . \ 
COREL ope boas MDRAOANMTMOOH DROaIWNMASTHOPE OOM 
OX OOOGORAK AKO ooon——— DOO OG wri ee ot ad od ed eed od mt ONIN 
DOODOCCOVOOC0C0000 OO Het ad ot 4 al 4 weeGaoctaaa od end owt 
DODDODOOODDODODOVDCDOCOCOOOCOOCCOOOCCCOCOCOO 


| 








| 





142 


— 
OY & 
>= 
~ 
ead —~ 
(ppifed ~~ 
— > »< 
Nt eS 
an —W 
| >*se -T 
i 2 3a 
4 = om <= ~ 
j nm LUA —_ 
| ~w TH —_ @ 
mA HH —— 
| N> we —S 
ae Ne Of <T 
mes Ne Io 
wN +#a j= 
Oe H> uu 
—Ht —t N= 
ll | < 
om ech 
i— Wet eww we 
— eer ee ~—w 
iN =~_J a> —_a 
> en ee _— om 
~N elu WN —_ ot 
Hommel, we if ww § ASW 
‘HumkmA Wt NN ee 
met Cl) mer Hh ~~ Li) 0 et ome te 
— Yam emu CO ww eww 
aet0O~ wa Zee Ne 
j-mt 3 Lh fw VO 
Wlmt mW Ei ee 
Ons #O H+ lvoe 
ow J) wer (0) om om ew —_—— om tf Ps 


IND e# ANY SAANIL a O@& 
tor eat GNA DeEUL II 

i PQENNSPOWOeHS Eee 
WW aww yy DN OW Ge 

Ut Seen eed coe | Mee? ST eel | ol 

[4 amet tant Grad teas oo tnd feet (CS tet tert tet om CN) OY” 
Zea IwOs NAINAw ere © 
meme YOdM=OOCDINAaWUZ 
NINN S&S I vOOOs>r- 00 rw 


um oO OO : 
N ST Ou 
Cl od “od 


DOANAPSTNOMOROANMTNO 
AMMAMMAMaAMMeTSITTS | 
ond G4 cond pad ad ed ge) ped el eed eed ed aed ed =) eed 


‘OOODDOOCOVOGCOCO0C0CO 


143 


TABLE V 


STATCR SCLLUTION 


F 
19) 
2 
7 
T 
0 





TABLE VI 


ROTCR SOLUTICN 


hr ODA Ht 
Ojbe eee @ 
IOONAKE 
PNMOS SU 
Aa ae et = 


WAOSN 
QDMOCWne 
Orato 
$~OTOOR 
eoeeee 


OOK aK 


crowurn 
AFournnr 
Tyroyreod 
m= OOM caoO 
woooedeai 
RJ] ee © @ 


ooco°o) 


PHONO O 
LONHOor 
Ioan rre 
mH ONOC 
Ly A tt HN) 
Bj @ @ ee @ 


oQooeco 


MOF 
WORN 
m— © @ @ @ @ 
OPrrinn 

ee Rene ute el ome 


PUNE O 
QM UWI 
uj @ eee e 
ear om 
bk PT TIN 
UN UIA LAU 


POOMDO 
NUNS ONO 
<—COmA 

eee e# 


COA 


DOAKAUD 
A DOWN 
CFWMown 

ge ade 


NOTTS 


NAIA 
HOMO 
ar or~r oO 
CN EAD EH 
eoecsc0ee 


CcCoOo0o 


MACAO 
Nato 
QOOxoS 
KDOWS 
LONWProO 
x eeee se 


CoOCCeH 


=m OON 
Osta On) 
GO) ee ee @ 
DPANUNDMOO 
mth A at 
MNF SU 


FRAC 


NPOITNS 
TFOPrNMNO 
be AUN WO ST et 
tli @ @ @ @ @ 
CORO wt ot ot 

OOM P= 


NO aa PO 
dU Onre 
LNG OMe 
a @©ee0ee 


I 


IANOONO | 
[on , 


NOR TOO 
aceonren 
NMAIMAIA 
O5eeeee 
be St st et 


Mp AIMO 
OxsFTuow 
“ON 0© ee @€ © @ 


ENA ONM | 


OM OA ER 


LAL LALA 


145 


UWI 
MH AMO+ 
“™ 0@© @ @© @ 
ek Onn Oh ole 9) 

TANIAIAIGC 

LA UA LALA 


OOTP 
WO OW t= 
eeeee 
NO Ofer 
bk AI OUCQUAION 
LA LAU LEA 


Vay oot tl) 
NMGMO ST 
e@e@ 8 e 
MOMMMN 
Spt ond ened ated ot 


| 


MFA COCUCN 
Pr ON 
he © e@ ee 
bm NAN SO 

MOMmMmMN 

LALA LALA 


prromo 
mMOROOO 


j= @ © @ © @ 


BOO 
od 


LOWUNaM 
MOST Oot 


ooreree 
Aanvnoeaor~ 


MN TIN 


WFHOO 
(a Al ood ph OL OX) 
aero 
LL PS UT OT ot ot 
CQOooccoo 


DOCO@ 


SAI OF end 
m= AIO UL\ ond 
bt ot ay fee) 
qevee 


hm POON 
LL GD 60 60 Cop 
eeee 


ddacc 


ud 
z 
= cath 


= 
Yn. 


FT-LB 


-3250-—-—- 14.663 


REF HORSE REF MOMENT 
POWER 


RSTAR 
0.40041 


Oe eR era eee me a oe mae nae meee em ee 


+s © came 


iG. 


INITIAL DISTRIBUTION LIST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library 
Naval Postgraduate School 
Monterey, California 93940 


Commander, Naval Air Systems Command 
Navy Department 
Washington, D.C. 20360 


Commander, Naval Ship Systems Command 
Navy Department 
Washington, D.C. 20360 


Capt. A. Bodnaruk, USN 

Naval Ship Systems Command (Code 6140) 
Navy Department 

Washington, D.C. 20360 


Office of Naval Research (Power Branch) 
Attn. Mr. J. K. Patton, Jr. 

Navy Department 

Washington, D.C. 20360 


Mr. R. Beichel 

Liquid Rocket Plant 
Aerojet-General Corporation 
Sacramento, California 95809 


Chairman, Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 93940 


Professor M. H. Vavra 
Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 93940 


Lt R. G,. Harrison, USN 

Naval Air Turbine Test Station 
1440 Parkway Ave. 

P.O. Box 1716 

Trenton, New Jersey 08628 


146 


No. Copies 


20 


SOCUnTy VleksSlUCaeuUon 


DOCUMENT CONTROL DATA - R&D 


(BSeeurity cleccitiecatiqn of tithe, bedy of abetract and indeming annotation muel be extored when the overall repart de alae se ied) 


1. ORIGINATING ACTIVITY (Cocpesate author) Ra REPORT SECURITY CLASSIFICATION 
3. REPORT TITLE . 
AN ANALYSIS OF SINGLE STAGE AXIAL-FLOW TURBINE PERFORMANCE USING THREE- 
DIMENSIONAL CALCULATING METHODS 



















4. OESCRI V& NOTES (Type of report and inclusive detec) 
Engineer Thesis, September 1967 
5. AUTHOR(S) (Laat name, fiset name, initial) 


HARRISON, Robert G. 


6. REPORT DATE Ja. FOTAL NO. OF PAGES 75. NO. OF REFS 
September 1967 144 9 


Ca. CONTRACT OR GRANT HO. Sa. CRIGINATOR*S REPORT NUMBER(S) 
















b. PAQJECT NO. 





10. AVAILABILITY /LIMITA TION NOTICES 


12. SPONSORING MILITARY ACTIVITY 
Naval Air Systems Command 
Navy Department 
Washington, D.C, 20360 


13. AGBTRAGT 


The method of turbine performance prediction developed by Vavra and 
Eckert has been refined in this analysis to realize more of the potential 
of the three-dimensional calculating methods, Mach number and rotor tip 
clearance effects on blade outlet angles and loss coefficient& have been 
localized rather than averaged over the blade height. An approximation 
for streamlane curvature has been used. 

Performance ‘curves were determined for two single stage axial-flow 


turbines located at the Propulsion Laboratory of the Naval Postgraduate 
School. Test results were available for one of the turbines. Agreement 
between predicted and experimental performance values was generally 
within 3 per cent. 





DD-1o2". 1473 UNCLASSIFIED 
1G7 Security Classification 


UNCLASSIFIED 


~ Security Classification 





KEY WOROS 


£NALYSIS 
SINGLE-STAGE 
kALAL@=FLOW 
TURBINE 
PERFORMANCE 
THREE=-DIMENSIONAL 
CALCULATING 
METHODS 





DD 72"",.1473 (sack) UNCLASSIFIED 


S/N 0101-807-6821 Security Classification A=31409 


148 























i il | lh Ty 


fi 





